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ABSTRACT 


Ionosphere-induced scintillation effects have been studied for 
many years. Only recently, however, have communications engineers begun 
to consider the detailed nature of the signal perturbation. Scintilla- 
tion activity was first characterized by a scintillation index and later 
by its first-order amplitude statistics. This gives sufficient infor- 
mation to determine the probability of occurrence of deep fades. 

However, if one considers the problem of optimizing receiver and 
signal design for combating scintillation effects, or of even evaluating 
them in a given system, it becomes clear that a much more detailed de- 
scription of ionospheric scintillation is necessary. Ideally, one would 
develop a channel model that describes the scintillation effects for an 
elementary waveform. Fourier techniques could then be used to analyze 
more complicated waveforms. 

The research reported herein provides a basis for developing such 
a channel model. The model would proceed from a parameterized functional 
description of the ionospheric spectral-density function. We have 
employed a power-law form with irregularity elongation along the direction 
of the geomagnetic field, but no attempt was made to determine the para- 
meter dependence on geophysical variables. 

We have characterised the diffracted field of a monochromatic plane 
wave by two complex correlation functions. For a gaussian complex field, 
these quantities suffice to completely define the statistics of the 
field. Thus, one can in principle calculate the statistics of any 
measurable quantity in terms of the model parameters. 
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The question of the validity of hypothesizing gaussian statistics 
for the complex field is an iinpox’tant one. At "best, gaussian statistics 
are a limiting form. Moreover, a considerable body of literature exists 
that questions the validity of this hypothesis. We therefore carefully 
analyzed the structure of the intensity statistics of several selected 
data sets for ionospheric and interplanetary scintillation. 

Our results show that the best data fits are achieved for intensity 
statistics derived under the gaussian-statistics hypothesis. The signal 
structure that achieves the best fit is nearly invariant with scintilla- 
tion level and irregularity source (ionosphere or solar wind) . It is 
characterized by the fact that more than 80 % of the scattered signal 
power is in phase quadrature with the undeviated or coherent signal 
component. Thus, the gaussian-statistics hypothesis is both convenient 
and accurate for channel-modeling work. 

To illustrate the technique of applying the channel model and the 
gaussian signal structure, we have computed the two-frequency correlation 
function for the complex field and its intensity as well as the wavelength 
dependence of the noimialized rms intensity scintillation index S^. The 
results are in agreement with reported UHP data, but they evidently 
cannot be extrapolated to L-band and higher frequencies on the basis of 
a simple power-law SDF. 

The technique of calculating phase statistics and xks nature of the 
limitation imposed by the implicit first Born approximation in our 
diffraction calculations are also discussed in this report. 
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I INTRODUCTION 


Ionosphere-induced scintillation effects are the main cause of signal 

degradation on satellite-to-earth radio transmissions at UHF and higher 

* 

frequencies. In the report hy Golden and IVolf (1973) , a need for a 
scintillation-effects channel model was clearly established. The research 
that is reported herein provides a basis for developing such a model. 

Channel modeling involves three tasks. Firstly, the ionospheric- 
ii'regularity structure must be determined. Secondly, a diffraction 
calculation must be performed to describe the electromagnetic-wave 
interaction with the scattering medium. Finally, the dependences of 
the model parameters on geophysical variables must be established. 

This work evolved from our efforts to extend and improve an existing 
scintillation model that we shall presently describe. Briefly, we have 

reformulated the diffraction theory so that it is more relevant to 

t'. 

systems analysis, and we have attempted to validate a gaussian signal- 
statistics hypothesis that greatly extends the usefulness of the theory. 

A , Background 

Scintillation activity is most commonly characterized by an ampli- 
tude (or intensity) scintillation index. Briggs and Parkin (1963) 

proposed four indices (S , S , S , and S ) and pointed out that the 

X ^ o ^ 

interrelations among them could be determined if the amplitude statistics 

were known. Their index is the standard deviation of intensity norma- 

4 

lized to the average intensity, and it is most convenient for theoretical 


References are listed at the end of the report. 
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computations. Briggs and Parkin (1963) performed a diffraction calcu- 
lation that relates to ionospheric and geometrical parameters. 

Bischoff and Chytil (I960) used the Nakagami distribution (Section II), 

-2 

whxch depends only on a parameter m = , to calculate Briggs and Parkin’s 

indices , and S^ as functions of S^. At the same time, because a 

large body of scintillation data existed on strip charts, AVhitney, Aarons, 
and Malik (1969) devised a scintillation index that could be determined 
from amplitude peaks. They showed empirically that their index was 
linearily related to S . 

Thus, a means existed for standardizing the various reported 
scintillation indices to S^, which could then be related to ionospheric 
parameters. In the Briggs and Parkin theory a random phase -changing 
screen with a gaussian autocorrelation function was postulated. The 
ionospheric parameters consist of an axial ratio, a transverse scale 
size, and an irnis phase variation. The rms phase variation is proportional 
to the ionospheric rms electron -density variation. 

In addition to these theoretical results, there existed a large 
body of published scintillation data. Hence, on 3 February 1971, NASA 
awarded a contract to SRI to develop a worldwide model for F-layer- 
produced scintillation (Fremouw and Rino, 1971, 1972). The model was 
based on a detailed enpirical formula for the ionospheric rms electron- 
density variation. Model constants were determined by using the Briggs 
and Parkin theory to fit carefully selected data sets. 

The scintillation model predicts an average scintillation index 
as a function of geomagnetic latitude, local time, season, and sunspot 
number. However, it was recognized that a scintillation index alone is 
insufficient for systems evaluation. At the very least, one needs the 
probability of occurrence of signal fades below a specified level 
(^Vhitney, Aarons, Allen, and Seemann, 1972). Hence a corollary study 
of amplitude statistics was conducted (Rino and Fremouw, 1973), 
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Our approach was to hypothesize gaussian statistics for the complex 
field and then calculate the amplitude probability "density function (PDF). 
The main difficulty lies with the fact that amplitude moments do not 
unic[uely determine the general gaussian PDF. The Rice, Rayleigh, and 
Nakagami PDFs are special cases that circumvent this problem. However, 
our initial calculations showed that none of these PDFs was accurate, 
although for some applications they would be usable. 

More disturbing was the fact that researchers in optics and radio 
astronomy had rejected the gaussian hypothesis altogether. They pre- 
sented both theory and data that showed the logarithm of amplitude to be 
normally distributed. ^Vhen gaussian statistics cannot be hypothesized, 
the channel -modeling problem is considerably more difficult. Hence, our 
first concern was to determine to what extent gaussian statistics were 
accurate. 

We also recognized a need for extending the diffraction theory 
so that the full temporal and spatial structure of the field could be 
determined. Ideally, the problem should be formulated within the frame- 
work of engineering systems analysis. This would make the results the 
kernel for the analysis of a broad spectrum of scintillation effects. 

On 11 November 1972, SRI received a second contract from NASA to 
carry out this work. The detailed objectives are summarized below. 

The results of the research are presented and discussed in Sections III 
through VI of this report, 

B. Objectives 

The objectives of this research were twofold. The first was to 
evaluate the gaussian hypothesis. The technique to be employed is 
described below. The second objective involved a purely theoretical 
effort also described below. The objectives were to be pursued simul- 
taneously with the hope that the results would be complementary . 
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Gaussian Signal-Statistics Hypothesis Test 


In an earlier document (Rino and Fremouw, 1973) we proposed a 
scheme for obtaining a family of general gaussian PDFs all with the same 
scintillation Index S^. The PDF family is parameterized by the ratio of 
the Fresnel-zone area to the square of the transverse scale size for a 
given effective axial ratio and incidence direction. Each PDF corresponds 
effectively to a different ionospheric rms electron-density value. 

Our first objective was to obtain a large body of data, compute 
histograms for each suitably stationary data segment, and obtain a best 
fit from the aforementioned PDF family for the loeasured S scintillation 
index. We then compared the fit to the corresponding log-normal PDF to 
determine which hypothesis (gaussian vs, log-normal) is more nearly 
correct, and to what extent. 

2, Computation of the Ionospheric Time-Variant Transfer Function 

In general, a system transfer function is the response (output) 
to a sinusoidal excitation (input). The ionosphere can be modeled with 
a suitable stochastic transfer function. Once the transfer function is 
determined, Fourier techniques can be used to determine the response 
to an arbitrary input. 

Thus, our second objective was to calculate the time-variant 
transfer function for a scattering layer that would provide a reasonably 
accurate representation of the ionosphere. The initial formulation was 
to be completely general and then specialized to a power-law spectral- 
densitv function. We proposed to make as few a-priori approximations 
as possible in order to test the various assumptions that are typically 
made. 

The theoretical calculations given define two complex correlation 
functions that suffice to completely specify the statistics of a complex 
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gaussiaix field^. Thus, in principle, one could calculate the statistics 
of any measurable quantity in terms of the model parameters. 

C. Organization of Report 

Because of the conflicting theories on signal statistics, we have 
devoted Section II of this report to briefly reviewing this subject, and 
to a discussion of the problem of testing the gaussian-statistics 
hypothesis from amplitude data alone. 

In Section III we describe in detail our data-reduction procedure 
and present the results of the data that we have analyzed to date. In 
Sections III-A and III-B we review the theoretical basis for analyzing 
histograms and applying goodness-of-f it tests. In Section III-C we 
review our method of obtaining a family of gaussian PD3?s each with the 
same first and second moments of intensity. 

In Section III-D we present the results of applying these procedures 
to ionospheric and interplanetary scintillation data. The results, 
which are discussed in Section III-E, support the conclusion that the 
gaussian-statistics hypothesis is more nearly correct than the log- 
noimal hypothesis. 

In Section IV we attack the general problem of channel modeling. 

In Section IV-A we discuss the philosophy of channel modeling and 
identify the quantities that must be computed. We then make a general 
computation for an arbitrary spectral-density fimction (Section IV-B) . 
This is followed by a specialization of the results for a power-law 
spectral-density function that admits a simple anisotropy (Section IV-C) . 

In Section IV-D we consider the effects of relative observer medium 
motion, which is the main source of the observed signal temporal 
fluctuations. In Section IV-E we discuss the method of evaluating the 
resulting iuteffi’als, and in the Section IV-F we characterize the general 
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properties of the transionosplieric channel and demonstrate them with 
some examples calculated hy numerically integrating the power-law 
formulas. 

Section V can he considered a synthesis of Sections IH and IV. 

In that section we apply the gaussian-statistics hypothesis to extend 
the complex correlations functions computed in Section IV to include 
intensity statistics. As examples we have computed the wavelength 
dependence of the S scintillation index and the two-frequency intensity 
correlation function. The results are in general agreement with published 
data. We have concluded the section with an outline of the method of 
making similar computations for phase statistics. 

In Section VI we have summarized the complete report and discussed 
its implications. We also consider the nature of the limitations imposed 
hy the first Born approximation and a method that we believe could be 
employed to extend oui’ results to allow for the effects of multiple 
scattering. 
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II 


SIGNAL STATISTICS 


Ratcliffe was aaiong tlie first researchers to suggest that radio- 
wave fading is caused by diffraction from ionospheric irregularities 
(Ratcliffe and JPawsey, 1933), By 1950, the theoretical basis for 
explaining the phenomenon was well developed (Booker, Ratcliffe, and 
Shinn, 1950; Briggs and Phillips, 1950; Hewish, 1952). The theory was 
restricted to single scattering, but Fejer (1953) had extended some 
results to include multiple-scattering effects. 

In the diffraction theories, the complex field is given as an 
integral over a random function that represents the scattering medium. 
Hence, by applying the central-limit theorem one can deduce limiting 
gaussian statistics for the complex field. However, evidently because 
of the influence of Rice's (1945) work, the additional assumption of 
random phase was usually made. This is well illustrated in Section 8.2 
of the popular survey paper by Ratcliffe (1956), 

Of course, other analytic techniques were exploited to analyze 
scintillation effects. Chandrasekhar (1952) used geometrical optics to 
calculate amplitude and phase scintillation for starlight. In 1960 
and 1961, however, English translations of the Russian monographs by 
Chernov (1960) and Tatarski (1961) became available. These were 
essentially complete treatments of wave propagation in random media. 
Their influence on the theoretical work that followed was considerable. 

V 

In the Russian works the Rytov solution to the wave equation was 

introduced. It predicts gaussian statistics for the logarithm of the 
\ 

comprex field rather than for the complex field itself. Moreover, it 
was originally thought to have a broader range of validity than the 



first Born approximation, which is implicit in the diffraction theories. 
However, this was later shown to he incorrect. What seems to he the current 
most widely accepted view is summarised in Barabanenkov et al. (1971), 

Our own results inject new evidence that favors the gaussian 
hypothesis over the log-normal hypothesis. Because this is a funda- 
mental question we have summarised the earlier arguments and the supporting 
evidence helow. 

A, The Log-Hormal Hypothesis 

The theoretical basis for the log-normal hypothesis is the Rytov 

solution to the wave equation. It gives the diffracted field E as 

®/e = where E. is the incident field and ilf is an integral over 

inc me 

a random function that represents the medium. Hence, an application of 

the central-limit theorem leads to the conclusion that log ^/e. has 

me 

gaussian statistics. 

We note, however, that when 1^[< < li ®/E. — 1 + ijr which is 

> me 

identical to the diffraction- theory solution. Hence, the Rytov solution 
appears to be more general. However, it has been more recently demon- 
strated that the condition |^|< < 1 is actually necessary for the validity 
of the Rytov solution (Brown, 1966), Nonetheless, the calculations 
require many approximations, and they are subject to some question. In 
view of this, theorists have considered the available data to help 
resolve the question. 

Most often cited are the optics data of Ochs and Lawrence (1969). 

Their plots of the logarithm of intensity versus cumulative gaussian 
probability show near linearity over three decades of intensity. More- 
over, their data fits for Rayleigh statistics are very poor by comparison. 
Nonetheless, their method of presenting the data tends to obscure small 
discrepancies , 
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More recently, Armstrong, Coles, and Eickett (1972) have analyzed 
radio-star intensity statistics for interplanetary scintillation data. 
Their data show clearly that log-noiraal statistics give subsT: ant tally 
better fits than do the corresponding Rician statistics. However, they 
do note statistically significant departures of their measured histograms 
from the corresponding log-no3cmal PDF. Thus, while their data are 
sufficient to reject the Rician-statistics hypothesis , they do not 
necessarily exclude the general gaussian-statistics hypothesis. 

In Barabanenkov et al. (1971) additional theoretical support for 
the log-noxmal hypothesis is given, however. By making a Markov 
approximation for the wave-medium interaction, it is possible to obtain 
an exact wave-equation solution for narrow-angle scattering. The 
calculated second moments evidently agree with the Rytov solution. 
Moreover, the calculated odd moments are finite. For gaussian statistics 
to be strictly applicable the odd moments must be zero. 

Even so, there have been no reported tests of the general gaussian 
hypothesis, for reasons discussed below. The difficulty lies with the 
fact that the gaussian hypothesis cannot be tested with intensity or 
amplitude data alone - 

B. The Gaussian-Statistics Hypothesis 

Our approach has been to hypothesize general gaussian statistics. 

The joint probability-density function of the real and imaginary parts 

of the complex field-say X and Y respectively— depends explicitly on 

2 

five parameters. These are the variances and means of X and Y, <J , 

2 ^ 
cr , 11 , "fl and their covariance C . None of these parameters can be 
y X y xy 

deteimined from intensity or amplitude data. However, in most experi- 
ments the undeviated component is effectively the phase reference, so 
that 1] = 0. 

y 
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2 2 

For Rayleigti statistics ct = ct , and C = Tj = 'll =0. Rician 

X y ' xy X y 

statistics permit a finite coherent component For these special 

PDFs, the moments of intensity uniquely determine the corresponding PDF. 

Hence, they are easily tested. The Nakagami PDF, which has been used 

extensively for ionospheric scintillation, has the same property. It 

-2 

depends on the parameter m = . Moreover, in some sense it is an 

approximation of the general gaussian amplitude PDF. 

However, as discussed in Rino and Fremouw (1973), the precise 
2 2 

relationship of m to a , a , C , and T] is not specified. \Vhen 

X y xy X 4 

is sufficiently small the Hakagami PDF is close to the Rician PDF. As 
approaches unity, the Nakagami PDF approaches the Rayleigh PDF. 

Between these extremes the Nakagami distribution differs substantially 
from the Rice distribution. Even so, we had no theoretical basis for 
expecting the Nakagami distribution to provide a better data fit than the 
log-normal PDF. 


In Rino and Fremouw (1973), we showed that by using a simple 

2 2 

diffraction-theory computation for c , cr , and C , and choosing T[ 

X y xy X 

2 2 2 1 | 2 , ^ ^ 
such that a + cr + 'll = { E. ; — that is, to ensure energy conser- 

X y X me 

vation — ^iiVe could obtain a family of gaussian PDFs all with the same first 

and second moments of intensity. Our initial results showed that the 

best data fit among the members of this family appeared to be as good 

as the log-normal PDF. The best-fit general gaussian PDF is distinguished 
2 2 

by the fact that cr > ^ or , thus explaining why the fits for Rayleigh 

V X 

^ i. 2, 2 \ 

and Rician statistics f'er - cr j are poor. Moreover, the condition 
22 ' ^ yf 

CT > > cr occurs when large-scale structure dominates the diffraction, 
y X 

Stated another way, we are observing a manifestation of the fact that 
the apparent scale size (squared) that produces amplitude scintillation 
is roughly equal to the Fresnel-zone area. 


10 



The idea of using general gaussian statistics to fit amplitude- 
scintillation data is not new. It is discussed in Wernik and Liszka 
(1969). Also, BiS‘’.hoff and Chytil (1969) discussed the possihle effects 
of on amplitude statistics. However, evidently the idea had not 
been previously tested, even though Bowhill’s (1961) calculations were 
directly applicable. 

In Section III we describe our data-reduction procedure and present 
our results, which show that general gaussian statistics as described 
give a data fit that is at least as good as, and generally better than, 
the correspohdlhg log-normal statistics. 
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III DATA REDUCTION AND ANALYSIS FOR SIGNAL STATISTICS 

In this section we shall make some general comments about histogx’ajns 
and goodness-of-f it tests. We shall then describe the theoretical basis 
for our da ta-r eduction pi'ocedure and present our results. 

The data that we have processed consist of digitized scintillation 
recoi’ds. The calibi-ated data samples constitute an instantaneous -power 
time sex'ies [P^^] for some pei'iod of time — say, T. The data are segmented 
into blocks of some conveniejit size for machine processing. For each 
data block an estimate of the average powei’ ^P) and the second moment 
of power (P ) is made. 

Generally, there are slow variations in the mean power level that 
can be compensated by renormalizing the data blocks to some convenient 
level. Once this is done, the data blocks are grouped into quasi- 
stationary periods. The stationarlty is subjectively evaluated by 

2 

obsei’ving changes In (p ) and carefully evaluating strip-chart recordings. 
It is Important for our analysis that genul.ne changes in the statistics 
do not get averaged out. 

Finally, we genei’ate histograms for each quasistationary data period 
and estimate the S scintillation index and, when enough data blocks ax’e 
available, its standard deviation. The histograms are recox’ded on digital 
tape foi’ coinpai’ison to various theoretical PDFs. Our main interest is in 
the best fit for the family of general gausslan PDFs that have the same 
first and second moments of intensity or power. 
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A. Histogram Computation 


Histograms are computed by dividing tbe power range into small 

intervals + ^^i^ ' counting the number of samples that fall 

into each interval. If the true PDF, f(P), is knom, one can calculate 

th 

the probability of one sample falling into the i interval. For 

small AP., Si f (p + AP./2) AP. . Prom P , one can calculate the 

X X 1 X 

histogram statistics. 


Indeed, for independent samples the number of samples that fall 
in the i interval (say, 0^) is governed by the binomial probability 
law 


P[0^ 


k] 



(III-l) 


It follows that 


(O.) = N P (III-2) 

X r. 

X 

and 

<0^> - <0.)^ = N P A - p \ (111*^3) 

" " ^x I 

Hence, the standard deviation of 0. is given as )/(N P_. ). Clearly 

X ^i ^i 

the uncertainty increases as P becomes small , so that one must make 

^i 

the usual compromise between resolution and precision. 

For convenience we have employed uniform increments. One would like 
at least five samples in the least significant interval to obtain good 
accuracy. However, this is not always possible for small N. A procedure 
for calculating confidence limits is given in Introduction to the Theory 
of Statistics by Mood (1950, pp. 233-235), when it is important to 
evaluate the significance of some segment of the histogram such as the 
tail. 

In our analysis, however, we have used a single measure of the over- 
all fit, as described below, 
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B, 


Goodness-of-Blt Tests 


Gooduess-of-flt tests are developed by devising a test statistic — 
that is, a measure of the discrepancy between the histogram (or some 
other quantity derived from the sample population) and the hypothesised 
probability law. Ideally, one would calculate the probability of the 
test statistic achieving a given level when the hypothesized probability 
law is correct. 


If the test is a good one, the test statistic will achieve a high 
confidence level only for sample populations that closely follow the 
hypothesized probability law. However, developing good test procedures 
is an art. Moreover, it is usually possible to calculate the limiting 
Gummulative distribution for the test statistics only as the number of 
samples approaches infinity. 

We have investigated the Kolmogoroff -Smirnoff (K-S) and the chi- 
square goodness-of-fit tests. The test statistics for the K-S test is 
the absolute value of the maximum departure of the estimated cummulative 
density function (CDK) from the hypothesized CDF. The limiting 
distribution for the test statistic has been determined. We found, 
however, that the K-S test is not sensitive to small parameter changes. 


The test statistic for the chi-square test is derived directly 
from the histogram. It is given explicitly as 



ClII-4) 


Hence, except for the normalization, it is a mean-square-error measure. 

- 2 ' ' 

Moreover, the test statistic X has a limiting chi-square distribution 
With! h degrees of freedom less one degree of freedom for each parameter 
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that must he estimated from the sample population (Bowker and Liberman, 
1950, pp. 458-461). 

We have found that Eq. (II 1-4) is very sensitive to changes in the 
gaussian parameters when they are varied as described below. Hence, we 
have used this test exclusively in our data analysis, 

C, Computation of Theoretical Gaussian PDFs 

The log-normal, Rayleigh, Rice, and Nakagami PDFs have the common 

2 

property that they are completely specified by the moments (P) and ) . 
By comparison, the general gaussian intensity PDF cannot be specified 
from intensity moments alone. One can show that for gaussian statistics, 


(P)= CT + T1 +11 

X y 


(III-5) 


and 


(P^> - <P)^ = 2/(<P> 


2 r r 

- CT ) 1 + ” 

a 


El “I 4 r lB| 

cos 2(^ - (zf) I + CT I 1 + 


1 4 r P 1 

) + CT 1 + — ;; — 

4 

-J L CT - 


where 


2 2 2 
o = CT + a 
X y 


(III-6) 

(III-7) 


1 1 ^ 

B exp {2i£l = a - CT + 2i C 
' X y xy 


ClII-8) 


and 


tan ^ = 


\ ■ 


(HI-9) 


2 2 'n <n 

We emphasise that none of the parameters cr , cr , C , H nor T] , can be 

X y xy X y 

determined from amplitude or intensity data alone. 


It follows from Eqs. (III-5) through (III-9) that there is an 
infinity of parameter combinations that give rise to the same moments {P) 
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/ 2 . 

and (P / . To make the problem oi testing the various possibilities 

tractable, we have used a diifraction-theory computation for B to derive 

equations for the necessary parameters. Suppose, for example, that 

2 

T] =0, and B is known. From Eq. (III-B) it follows that S can be 
y 4 

written as 


2 2 4 

= 2 o g + cr (g - 2g ) 
4 1 2 1 


(III-IO) 


where 


/1b|\ 

= (~) 

\ <^ / 


cos 2 Q 


(III-ll) 


and 


-1+ ^ 


(III-12) 


2 2 

Now, if we measure S , o' can be determined from the formula 


cr 


= <P) 


r 2 2-]^/^ 

"■ L®1 4j 


[gg - 2g^] 


(III-13) 


which gives the positive root of Eq. (III-IO). Once cr is determined, 
the gaussian parameters follow fromEqs. (111-5), (III-7), and (III-8) 
as - 



2 

= [<P> - o-^]' 


(111-14) 


2 

CT 

X 



1 + 


Re{B} 

2 

C 


(m-15) 
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1 2 
2 


(III-16) 


1 2 

2 ^ 2 
CT 


(III-17) 


A simplified diffraction calculation with a gaussian autocorrelation 
function for the irregularity structure gives the result 

2 f/ 2 \^ / \2 

E = - cf fl“tan tan Sec 9j + ff^ tan tan U^j 


/ f tan U + f„tan U 
1 - 1/1 12 2 
X exp - 1 tan I — 

\ 1-tan U tan U sec 9 
_ \ 1 2 


(III-18) 


tan = — 

P 


CIII-19) 


tan U„ = Z 
2 


Cin-20) 


2 2 

f = 1 + tan 0 cos cp 
1 


(III-21) 


2 2 
f = 1 + tan 0 sin 
2 


cm- 22 ) 


The parameter Z is the ratio of the Fresnel-zone area 


sec 9 j 


to the 


transverse scale size squared. The parameter 6 is related to the axial 


ratio, a, and the geomagnetic dip angle, ijf , hy the relation 
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3 = 


( 2 

[a cos 


i]f + 


2 

sin ijf 


) 


1/2 


cm-23) 


The angles 0 and cp are the incidence angles at the ionospheric penetration 
point (see Rino and Fremouw, 1973) , 

From Eq. (III-18) we see that as Z varies from zero to infinity, 

2 

B varies from -a to zero. Hence, from Eqs. (III-IS) and (III-16), we 

2t 2 

see that for small Z, cr > > cr . Only in the limit as Z approaches 

^ 2 2 2 \ 

infinity are Rician statistics — cr /2.'j achieved. Hence, for 

a given value of S^, we can change the form of the intensity PDF 
considerably by varying the single parameter Z. We must, of course, 
first fix the axial ratio a and determine the angles 0, tp> and ilr. 

Examples of two families of curves obtained in this manner are 
shown in Figures Ill-l(a) and Ill-l(b), For comparison we also show the 
corresponding Nakagami and log-normal PDFs, One can see clearly that 
there are significant differences for small Z values, particularly for 
the larger scintillation-index curves. Hence, if the data show similar 
differences, they should be detected by careful processing. The 
question of the validity of Eq. (III-IS) will he deferred until the 
results are presented. 


D. Data Analysis 

1. NASA OAQ-2 Data 

Simultaneous data at 136 and 400 MHz are transmitted by the 
OAO-2 satellite, which is in polar orbit. The purpose of the trans- 
mission is to telemeter data on the 400-MHz channel. The satellite was 
routinely observed at the NASA tracking at Quito, Ecuador during 1969 
and 1970, Data were selected for scintillation analysis when operator 
logs indicated that difficulties were being encountered in applying 
normal operating procedures. 
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LA-2273-12 


FIGURE m-1 FAMILY OF GAUSSIAN PDFs PARAMETERIZED BY Z COMPARED WITH 
THE CORRESPONDING NAKAGAMI AND LOG-NORMAL PDFs. For large 
Z, the gaussian PDF is very nearly Rician. 
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Data for eight representative passes , selected as described 
above, have been processed. Four data channels (nUJUbered 1, 3, 5, and 
7) , consisting of horizontal and vertical polarizations at 136 and 400 MHz 
respectively, were recorded for each pass. The data base is summarized 
in Table III-l. The data were generally recorded at low elevation angles, 
and they were characterized by very intense scintillation. Indeed, the 
S index was often greater than unity at 400 MHz, 


Table III-l 


NASA OAO-2 DIGITAL-DATA TAPES 


■ 

Analog 

Number 

Date 

GORE Tape Numbers 

Time (UT) 

1 


Channels 
3-7, 5-7 

F 

a 

1 

1 

E 

1040 

1082 

10/1/69 

10/14/69 

BW7975 

BW9421 

BW9417 

BW9433 

0951-1006 

0314-0329 

q 

u 

i 

1085 

10/15/69 

BW9419 

BT4556 

0235-0257 

n 

0 

X 

1092 

10/17/69 

BT4578 

BT4587 

0304-0319 

S 

P 

1489 

03/3/70 

BW9544 

BX0487 

0055-0110 

V 

i 

li 

g 

15S6 

03/28/70 

BX0481 

BWSS90 

0220-0233 

E 

q 

u 

1660 

03/30/70 

BV9545 


0233-0248 

i 

n 

1665 

03/31/70 

BX0239 

BX0240 

0340-0356 

X 
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The measured S scintillation indices for the fall equinox 
4 

data are summarised in Table III-2, Each data block encompasses approxi- 
mately 15 seconds of data. A most disturbing fact is that the statistics 
are dissimilar at the same frequency for different polarisations. The 
single exception is the 1 October 1969, 400“MHz data for Blocks 3, 9, 
and 10. 


Table I I 1-2 

VALUES FOR OAO-2 FALL EQUINOX DATA 
4 


(a) 1 October 1969, 0959:15 UT 


Block No. 

BW 7975 File 1 

BW 9417 File 2 

136 MHz Horiz. 

136 ftttiz Vert. 

400 MHz Horiz. 

400 MHz Vert. 

1 

1.1789 

2 . 1233 

0.4071 

0.6052 

2 

1.2252 

2.7863 

0,7741 

2.6971 

3 

1.0005 

2 . 4049 

0.6281 

0.2355 

4 

0.9925 

1.8213 

0.7195 

0.5011 

5 

1.1087 

2.0123 

0.1931 

0.2767 

6 

1,3726 

2,2537 

0.4766 

0.3655 

7 

1.05124 

3.2596 

0.5740 

0.6639 

8 

1.5115 

4.7579 

0.9410 

0.9278 

9 

1 . 5744 

3.1266 

0.7641 

0.7553 

10 

0.7745 

1.4798 

0.2699 

0.2757 
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Table III-2 (contirLued) 


(b) 14 October 1969, 0316:04 UT 


Block No. 

BW 9421 File 1 

BW 9433 File 2 

136 MHz Horiz. 

136 MHz Vert. 

400 MHz Horiz. 

400 MHz Vert. 

1 

0.9339 

0.8759 

0.8817 

1.0477 

2 

1.9755 

2 . 5796 

0.7025 

0,8122 

3 

1.S632 

2 .1910 

0.8834 

0.6527 

4 

2.7457 

3.8259 

0.9971 

0.8230 

5 

1.3776 

1.9358 

0.9317 

0.8113 

6 

2.4834 

2 . 9448 

1.9673 

1 . 7433 

7 

2.2778 

2.8537 

1.5217 

1.2464 

8 

1.9309 

2 .2038 

1.0278 

2.4421 

9 

2.0347 

2.3987 

3.6043 

1.3116 

10 

1.7913 

1.8859 

1.5389 

0.7388 

11 

2.7145 

2.9625 

0,7257 

0.5043 

12 

1.9286 

2.1612 

0.5467 

0.4050 

13 

1.8836 

2.4187 

0.2795 

0.2372 

14 

1 . 3492 

1.0591 

0.6565 

0.5001 
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Table III-2 (continued) 


(c) 15 October 1969, 0244:03 UT 


Block No. 

BT 9419 Pile 1 

BT 4556 File 2 

136 MHz Horiz. 

136 MHz Vert. 

400 MHz Horiz. 

400 MHz Vert, 

1 

2.8125 

4,4357 

1.0061 

0.5967 

2 

1.1833 

2.7865 

1.0472 

0.6532 

3 

2.2734 

2 . 8146 

0.7737 

0.8753 

4 

2.0565 

3.1683 

1.5404 

0.8502 

5 

1.9092 

2.7466 

1.0703 

1.3409 

6 

1.8667 

2,7279 

1.2106 

0.8164 

7 

1.9109 

2.2961 

0.8316 

0.7684 

8 

1.6300 

1.5960 

0.5620 

0.7638 

9 

1.665S 

1.5764 

1.0015 

1.4876 

10 

1.7062 

1.8075 

0.3899 

0.4950 

11 

2.0086 

1.9361 

0.2221 

0.3846 

12 

1.5782 

1.4455 

0.2316 

0.3820 

13 

1.5835 

1.4111 

0.3557 

0.4869 

14 

1.9054 

1.5007 

0.3103 

0.4052 

15 

1.9456 

2.3084 

0.3961 

0.5007 

16 

2.1826 

2.0509 

0.5640 

0.5809 

17 

1.8364 

1.8634 

1,0538 

1.3310 
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Table III-2 ^concluded) 


(d) 17 October 1969, 0308:01 UT 


Block No. 

BT 4578 File 1 

BT 4587 File 2 


136 llHz Horiz. 

136 MHz Vert. 

400 MHz Horiz. 

400 MHz Vert, 

1 

1,3234 

1.1452 

0,3656 

0.3277 

2 

1,8886 

1.4484 

0.2815 

0.3012 

3 

1,5665 

1.2048 

0.4251 

0.4375 

4 

.1.8715 

1.2709 

0.3925 

0.4126 

5 

2.7674 

2 . 8522 

0.5915 

0.7938 

6 

7.4664 

5.0305 

1,5980 

1,0067 

7 

5.6374 

4.3200 

0 . 5403 

1.7873 

8 

3.7195 

2.9596 

0.3623 

0.5829 

9 

3,7616 

3.0175 

0.4199 

0,5931 

10 

2 . 8402 

2.9652 

0.3279 

0.3569 

11 

1.5799 

2.1713 

0.2110 

0.0982 

12 

8.4029 

1.2729 

0.5618 

0.1661 

13 

6,3948 

17.7023 

0.3794 

0,3641 

14 

3 . 5439 

3.9836 

0.9426 

1.0086 

15 

6.3153 

1 . 0640 

0.5162 

0.4097 

16 

1,0136 

1.6912 

1.2929 

0.8343 

17 

2.0756 

2.3189 

0.7393 

0.7099 

IS 

2 . 9416 

3.4522 

1.7486 

0.9880 

19 

2 . 9619 

3,9994 

1.5129 

0.9516 

20 

6.5558 

7.9449 

3.1781 

2,9251 

21 

8.0384 

1.8612 

1.3810 

1.5783 

22 

13.5004 

7.7942 

1. 0309 

1.0226 

23 

9.3590 

5.6452 

1.4560 

1.1452 

24 

2.7745 

3.2785 

1.8369 

1.5030 

25 

1.6137 

1.7660 

0.2078 

0.2250 

26 

1.9871 

1.8546 

0.7048 

0.4457 

27 

2.2410 

1,7230 

0.4764 

0.5736 

28 

4.9856 

2.5434 

0.7316 

0.9798 

29 

5.6505 

3.0766 

0.8343 

1.2218 

30 

7.3730 

3.. 2692 

0.6834 

0.9621 












NUMBER OF COUNTS (TOTAL 6000) NUMBER OF COUNTS (TOTAL 6000) 


The a'oove mentioned fact, coupled witli the generally strong 
scattering conditions and the short stationary periods due to satellite 
motion, renders these data unusable for pur purposes. We note that a 
gaussian field cannot produce an index greater than tjz. Nonetheless, 
we have processed one data set to illustrate the difficulties. 

For the 16 October 1969 data, only the first of the pair of 
digital tapes that constitute one pass was available. Hence, only the 
136-MHz data were complete. Nonetheless, they display generally weaker 
scattering than for the data summarized in Table III-2, The processed 
data for the horizontally polarized antenna is shown in Figure I I 1-2. 



°o 1.0 2.0 3.0 4.0 

NORMALIZED INTENSITY 



FIGURE a [-2 LOG-NORMAL AND BEST-FIT GAUSSIAN PDFs FOR NASA 136-MHz OAO-2 
DATA — ANTENNAS HORIZONTALLY POLARIZED 
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NUMBER OF COUNTS {TOTAL BOOO) NUMBER OF COUNTS (TOTAL 6000) 


The best-fit gaussian PDF is clearly a better fit to the data than the 
corresponding log-normal PDF, However, for the corresponding vertically 
polarized data shoivn in Figure II 1-3, neither PDF provides a good fit. 




LA-2273-2 


FIGURE 111-3 LOG-NORMAL AND BEST-FIT GAUSSIAN PDFs FOR NASA 136-MHz OAO-2 
DATA — ANTENNAS VERTICALLY POLARIZED 

Because such behavior cannot be reconciled with any currently 
available weak or strong scatter theory, we did not pursue these data 
further. Rather, we have concentrated on synchronous-satellite data 
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taken tinder conditions of weaker scattering. However, tke QAO-2 data 
are interesting in their own right. They have "been analyzed (Blank 
and Golden, 1973) and the polarization effect noted. 

2. HAS A ATS “5 Data 

Six data sets for the synchronous satellite ATS-5 recorded at 

(i 

Lima, Peru during Hoyemher 1971 have been processed. Among the six data 
sets only one produced processable data. The data tapes and the 

I 

difficulties encountered are summarized in Table III-3 below. 


Table 1 1 1-3 
HASA ATS-5 DATA 


. GORP Ho. 

Date 

Comments 

B07445 

11/17/71 

Short record 

B05334 

11/02/71 

Bad tape 

B07830 

11/11/71 

Good data 

B07805. 

11/12/71 

Bad calibration (possibly recoverable) 

B06819 

11/08/71 

Wrong format 

B063 }7 

11/04/71 

Wrong format 


Two channels of data were recorded with spaced receivers 
separated by 1200 ft. along an east-west baseline. The receivers were 
located at Lima, Peru. These data are identical in format to those 
presented by Hino and Fremouw (1973) . The two data channels should have 
nearly identical statistics. The data tapes B06S19 and B06327 had very 
different data in the two recorded channels. Hence, we did not process 
these data further. 

One complication with the ATS-5 data lies with the fact that 
the satellite is spin- stabilized. The antennas are not completely 
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despun, so that the data show a slight periodic modulation at approximately 
0,16 Hz. The effect is to smear the histograms, since data from different 
parts 01 each cycle have different means. We did not attempt to correct 
the data fcr chc satellite spin because of the increased complexity and 
cost. 

In Figure III-4 we show the histograms for the 11 November 1971 
data together with the corresponding log-normal PDFs. The measured S 
indices Were 0.37 and 0,34 for Channels 1 and 2, respectively. There 
is clearly a statistically significant difference between the measured 
histograms for the two channels. The near 10% difference in the 

indices is larger than their individual standard deviations of ~ 2%. 

2 

Moreover, the respective X parameters differ by more than a factor of 4. 

Evidently the two receiver characteristics are different. For 
example, the satellite spin is slightly more pronounced in the Channel-2 
data. This was inferred by noting that the relative intensities of the 
spectral lines due to the satellite spin* are approximately 1 dB higher 
for the Chamiel-2 data. 

The results of applying the curve-fitting procedure that we 
have described reflects the differences between the two channels. In 
Figures III-5 and III-6 we show the best-fit gaussian PDFs, For the 
Channel-1 data the improvement is considerable. We see that the "rise 
time" for the log-normal PDF is too slow. Moreover, it underestimates 
the peak of the histogram. The best-fit gaussian PDF completely corrects 
these discrepancies (cf . Figure III-4), The best-fit parameters are 
summarized in Table III-4, We see that for both channels, less than 20% 
of the scattered power is in-phase with the undeviated competent. 


See Figure II 1-7. 
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NUMBER OF COUNTS 


DATE = 11/11/71 TAPE = B07830 FILE = 1 



FIGURE !ll"5 BEST-FIT GAUSSIAN PDF FOR NASA 134-MHz ATS-5 DATA — CHANNEL 1 
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Table I I 1-4 


PARAMETERS FOR MSA ATS -5 DATA 
(3 = 5.0, 0 = 31.04°, 05 = 110.8°) 



Channel 1 

Channel 2 

s 

0.3724 

0.3421 

4 

■ 


2 

X 

719.0 

816.0 

Z 

0.590 

0.880 

2 

a 

0.2162 

0.1620 

2=1= 



O' 

0.0811 

0.1357 

X 

« 


|b|* 

0.8892 

0.8082 

Z.B 

19.00° 

25.60° 


♦ 

Normalized to incoherent component 
rms power. 


For the Channel -2 data the log-noimal PDF remained a slightly 
better fit than the best-fit gaussian PDF. Hence, the data are not 
mutually consistent. However, as noted above, the more pronounced spin 
effect may be cause for rejecting the Channel -2 data. 

To complete analysis of the ATS-5 Lima data we have computed 
the power spectrum of the intensity fluctuations. The results are shown 
in Figure III-7 together with a family of power-law SDFs [sq. (IV-27)] 
for different spectral indices y. The y = 1.5 curve, which is evidently 
the best fit, is close to the y = 1.45 value inferred from satellite in- 
situ measurements (see Section V-D) , However, the receiver noise tends 
to force this result, and a more careful analysis should be performed. 
The arrows point to the satellite spin lines. 
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RELATIVE POWER 



FIGURE III-7 INTENSITY POWER SPECTRUM FOR NASA ATS-5. Arrows point to satellite spin lines. Solid curves are for 
power-law analytic forms. 



3. AFGHL ATS-5 Data 

The Air Foi’ce Caiubi’idge Research Laboratory has processed a 
large volume oi scintillation data for amplitude statistics. Their 
primary interest was in characterizing the average statistics for 
relatively long periods of ti:iie. They found the Nakaga:!ii distribution 
to be well suited for that particular task (^Vhitney, Aarons, Allen, and 
Seemann, 1972). 

Because of our mutual interest in intensity statistics, 
arrangeJiients were made for SRI to obtain and process some representative 
samples of the AFCRL data. Our analysis of the data was complementary, 
in that we are interested in a precise computation of the statistics for 
short periods. In this subsection we svimmariae the results of processing 
sixteen hours of data recorded at Hamilton, Massachusetts during April 
and May of 1973. Both VHF and UHF signals were recorded, although the 
UHF scintillation was generally too low to be reliably processed. 

The data are grouped into five quasistationary segments, as 
described in the introduction to this section. The data base is summarized 
in Table I I 1-5 , Data Sets 2 and 4 contain only 1000 samples each. Hence, 
their accuracy is generally poorer than for Data Sets 1, 3, and 5. For 
each data set we have obtained the best-fit gaussian and the corresponding 
log-normal PDFs, 

The parameters fox* the VHF (136 Mia) data are summarized in 

2 

Table I I 1-6, We first note that for each data set the gaussian X value 
2 

denoted as X achieves a lower value than the corresponding log- 

2 gauss g 

normal X value denoted as X • Moreover, the difference is lai’gest 

LH 

for Data Sets 1, 3, and 5, which are statistically most significant. 
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Table III-5 


SUMMARY OF AFCRL ATS -5 DATA 


Data Set 

Date 

No . of Hours 

Sample Rate 
(Hz) 

S (VHF) 
4 

1 

4/12/73 

2 

0.14 

0.5527 

2 

4/12/73 

1 

0.14 

0.3442 

3 

4/15/73 

3 

0.14 

0,2664 

4 

4/15/73 

1 

0.14 

0.1356 

5 

5/18/73 

4 



0.31 

0.2829 


Table I I 1-6 

PAMHETERS FOR VHF DATA 
(P = 3.2, 0 = 55.20°, ^ = 40.21°) 


Data 

Set 

^4 

2 

2 

X 

gauss 

2 

a 

2 

CT 

X 

|b| 


Z 

1 

0.5527 

175.4 

160.5 

0.3501 

0.0971 

0.8780 

23.4° 

mm 

2 

0.3442 

113.9 

113.8 

0.1494 

0.1018 

0.7879 

31.1° 


3 

0.2664 

86.5 

56,67 

0.1091 

0.1307 

0.8329 

27.5° 


4 

0.1356 

82.7 

75.64 

0.0145 

0.3171 

0.5481 

48.1° 


5 

0.2829 

169.5 

107.2 

0.1362 

0.1005 

0.8734 

23.8° 

0.41 
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We also note that for Data Sets 1, 3, S, the Z parameter lies 

2 

between 0.4 and 0,5. The cori'esponding cr parameter’s show that less than 

X 

20% of the scattered power Is in phase with the steady or undevlated 
signal components. The equiprobahillty ellipses for the joint PDF of 
X and Y are chax'acteriaed by |b| and /iB. Indeed, one can show that the 

1 + [b1 

ratio of major to minor axes is given as — ^ — — r, and the orientation 

1 - |b| 

angle is 1/2 <iB. 

Hence, the best«fit gaussian PDF has an axial ratio greater’ 
than 10 (not to be confused with the iri’egularity axial ratio a) and it 

Q 

is within 10 of being perpendicular to the undeviated signal component. 

The strength of the scattered component is the main variable. We note 

2 

that there is no obvious breakdown of the theory for cr as large as 0.35. 

It is worth noting that while we have discounted the a'i.O-B data, the PDF 

2 

in Figure I I 1-2 (a) corresponds to a o parameter larger than 0.9, 

In Figures III-8, III-9, and III-IO we show the log-normal PDF 
and the best-fit gaussian PDF for Data Sets 1, 3, and 5. We see in each 
case that the log-noimial statistics achieve a slower rise time, and then 
underestimate the peak of the histogram. This pattern was previously 
observed in the NASA ATS-5 data at least for the Channel-1 data [cf . 
Figures III-4(a) and III-5(a)]. 


In Table 1 1 1-7 we have sunmiarlzed the UHF parameters for the 
three data sets that showed the strongest scintillation. The parajiieter 
Tj is the spectral index defined in Section V [Eq. (V-3)]. We first note 


2 ^2 VHP 

that theory predicts e “ A. . However, for Data Set 1, -r — = 11.9, 


VHF 


UHF 


whereas = 9.18, We believe the discrepancy to be due to noise 

UHF 

contajiiination of the VHF data. The spectral-index value is also somewhat 
smallei’ than we expect (see Section V-A-1) . 
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TAPE = 7045 FILE “ 1 




FIGURE m-8 BEST-FIT GAUSSIAN AND LOG-NORMAL PDFs FOR AFCRL DATA — SET 1 
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Table I I 1-7 

PARAMETERS FOR UHP DATA 
(3 = 2.1, 6 = 31,04°, 0 = 10.8°) 


Data 

Set 

^4 

2 

2 

X 

. gauss 

2 

cr 

2 

CT 

X 


Lb 

Z 

T] 

1 

0.0831 

123,8 

123.8 

0.0293 

0.0464 

0.9431 

15. 9° 

0.25 

1.76 

2 

0.0714 

61.87 

61.87 

0.0126 

0.0971 

0.87S0 

28.4° 

0.4 

1.46 

5 

0.0634 

SO.OQ 

65.61 

0 . 0023 

0.43 

0.32 

65,4° 

3.2 

1.39 


Ttie diserepaucies become larger as the UHF index becomes 
smaller. Moreover, we observe that the best-fit Z parameter steadily 
increases. That is, the statistics are tending to Rician as we would 
expect if the noise contribution were dominating the statistics. We do 

note, however, that for each data set the best-fit gaussian PDF achieves 

3 2 

a lower X value than the corresponding log-normal X value. 

4 . UCSD Interplanetary Scintillation Data 

As a final example we present some interplanetary radio-star 
scintillation data talien at the University of California Solar Observatory 
located near San Diego, California. The observations were made at 
75 MHk. The UCSD data processing for amplitude statistics is discussed 
by Armstrong, Coles, and Rickett (1972), which we referred to in 
Section II-A. 

Arrangements were made to apply our data-r eduction procedure 
to some representative UCSD data. We assumed normal incidence and 
iostropic irregularities. In the data reduction we varied the transverse 
scale size rather than Z directly. The wavelength is known, and the mean 
distance between the earth and the sun was used for the distance para- 
meter z. 
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Three typical results are shown in Figures Ill-llCa), (h) , 
and (c) . In all cases the hest-f it gaussian PDF fit better than the 
corresponding log-normal PDF. The transverse-scale parameter varied 
between SOO and 1000, which is comparable to the square root of the 
Fresnel-aone area for 75 ^'[Hz at the distance of the sun. Thus, the best- 
fit S-parameter value is near or less than unity just as it was for the 
ionospheric scintillation data. Indeed, the strihing feature in these 
data is the similarity of the histograms to those obtained for the 
ionospheric scintillation data. We observe again that the discrepancy 
between the data and the log-normal PDF is systematic in that the log- 
normal PDF achieves too small a rise time, and then underestimates the 
peak of the histogram. 

The similarities between the ionospheric and interplanetary 
scintillation strongly si’gge^t that certain features of the scattering 
are universal. This posa-ioility is discussed in the next subsection, 

E, Summary and Conclusions 

In this section we have described a procedure for obtaining a family 
of probability-density functions for intensity based on the assumption 
of gaussian statistics. Each member of the family has the same first and 
second moments of intensity. By comparison, the log-normal and Nakagami 
PDFs are uniquely determined by the first and second moments of intensity. 

In Section III-D we applied the procedure to obtain a best fit to 
measured histograms from ionospheric and interplanetary scintillation 
data. The results support the general conclusion that the best data 
fit among the members of the gaussian family is significantly better 
than the corresponding log-normal PDF. The implication is that the 
diffracted field rather than its logarithm has gaussian statistics. 

The hest-fit gaussian PDFs display a considerable degree of 
similarity. We find that independent of the value of S , the magnitude 
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PROBABILITY (arbitrary units) PROBABILITY larbitrary units! PROBABILITY [orhitrary units) 




0 1 2 3 A 

NORMALIZED INTENSITY 


LA-2273-17 


FIGURE 111-11 BEST-FIT GAUSSIAN AND LOG-NORMAL PDFs FOR UCSD 
INTERPLANETARY SCINTILLATION DATA 


o 

of B is generally larger than 0,8, while its angle is less than 30 . 

This implies that moi-e than 80% of the scattered power is in phase 
quadrature With the imdeviated signal component, and is finite. 

These results readily explain why Rician statistics generally 
provide a poor fit for scintillation data. The hest-fit gaussian PDF in 
general tends to he more peaked about its mean than either the corresponding 
Rice, Rakagami, or log-normal PDFs. As a ‘practical engineering matter, 
the latter PDFs will give more conservative probabilities for a given 
index. As an example, we have computed fade margins for a given 
index and different PDFs by setting B = 0.85 which seems to a 

good median value* The results are shown in Figure III-12 



LA-2273-18 

FIGURE 111-12 FADE-MARGiN COMPARISON FOR GAUSSIAN. LOG-NORWIAL, 
AND NAKAGAMl PDFs 
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Finally, we consider tlie validity of EQ. We know 

from the shape of the intensity power spectrum [see Figure III-71 that 

a gaussian autocorrelation function is incorrect. Indeed, if one uses 

e point as a measure of spectral width, the gaussian PDF corresponds 

2 

y 03 limit. To see this, consxder that e = Ixm 

What we have observed suggests that any diffraction model that 
gives rise to the appropriate B value for some set of parameters will 
achieve a near-optimum data fit. It seems unlikely that one could 
obtain sufficiently accurate data to resolve fine details even if a 
unique solution did exist. 

We shall show in Section IV that independent of the particular 

spectral-density-function shape, one can always achieve the condition 
2 2 

or > > cr for sufficiently small Z. The transverse scale parameter 
y X 

can be interpreted as an "outer irregularity scale" following the 
terminology of turbulence theory. A small value of Z will always be 
achieved if the outer irregularity scale is sufficiently large. Indeed, 
for incompressible-fluid turbulence, the outer scale is infinite. Our 
results are consistent with the assumption of a similar result for 
magnet ohydrodsmamic fluids. 

In closing this section we note that one possible way to 
reconcile with theory the result that the field statistics are gav.ssian 
rather than log-normal (see Section II) is to reject the Markov assumption. 
It has been shown that the Markov assumption is sufficient for the 
validity of the Rytov solution, Tatarski (1971) has discussed the 
validity of the Markov approximation and sho^vn that it is consistent 
with iostropic irregularities, but not necessarily with anisotropic- 
irregularity strictures. Since the latter is the case for ionospheric 
irregularities we have some basis for believing that the Markov assumption 
is in fact too stringent for radiowave scintillation. 
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IV THEORETICAL BASIS FOR CHANNEL MODELING 

In the introduction to this report we briefly discussed the general 
problem of channel modeling. In this section we shall present in detail 
theoretical calculations that form a basis for a scintillation-effects 
channel model. 


A. Introduction 

A functional block diagram of an artificial earth satellite commu- 
nication link is shown in Figure IV-1. For our purposes, the modulation 
is a complex voltage applied for some period of time— say, T milliseconds. 
Demodulation is the process of removing the complex envelope from the 
received signal. If the transionospheric channel were perfect, the 
demodulated signal would be an exact replica of the modulating waveform. 




LA-227 3-5 


FIGURE IV-1 FUNCTIONAL BLOCK DIAGRAM OF SATELLITE-EARTH COMMUNICATION 
SYSTEM 


We are not concerned here with additive noise or interference, 
since even in their absence the demodulated signal is not a perfect 
replica of the modulating signal. A wave packet suffers temporal 
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dispersion owing to the fact that different frequency components propa- 
gate with slightly different velocities. Moreover, raypaths are curved 
so that a "pencil” of rays suffers spatial dispersion as well. 

There are several other effects, such as Faraday rotation and 
absorption, to name two. These effects are characterized by the fact 
that they can, in principle, be computed from the deterministic component 
of the ionospheric index of refraction, which depends explicitly on 
frequency and position (Lawrence, Little, and Chivers, 1964). 

Our interest here is in the effects produced by the irregular 
departures of the ionospheric index of refraction from its mean value. 

Thus, we must make a somewhat arbitrary division between the slowly 
varying (spatially as well as temporally) component, which we shall 
assume to be deterministic, and the residual, which we assume to be 
representable by a suitable random function. 

For simplicity, we shall ignore the deterministic effects, so that 
in the absence of index-of-refraction irregularities, the signal free- 
space-propagates from the transmitting antenna to the receiving antenna . 

There is no serious loss of generality in neglecting the deterministic 
effects; some can be eliminated by operating procedures. For example, 
Faraday-rotation effects can be minimized by using circularly polarized 
antennas. 

For engineering applications we seek a description of the terminal 
behavior of the system shown in Figure IV-1. Indeed, a fairly well developed 
theory of signaling in fading dispersive channels exists. One of the most 
basic measurements is a determination of the system frequency response. 

For the transionospherlc channel, a monochromatic signal emerges with 
a stochastic component that is slowly varying compared to the frequency 
of the incident wave. 
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Tliere are two quantities to be determined — namely, the fraction 
of the incident enex'gy that is randomized, and the statistical structure 
of the random component. Our analysis will be I'estricted to weak 
scattering, and we shall first calculate the complex second moments 
of the field. In Section V we shall apply the properties of gaussian 
fields to calculate intensity statistics. 

Following oux’ "blaok-box” approach to channel modeling, we would 
calculate the complex second moments of the demodulated voltage 
vCt) = ^v(t)) + 6vCt).. The technique has been discussed in detail 
by Bello (1963). We represent the response of the system shown in 
Figure IV-1 to a sinusoid of frequency f, as h(t;f), a stochastic 
function of time t and frequency. Then, the response to an arbitrary 
signal whose specti’um is x(f) can be calculated from the integral 

v(t) = J* x(f)h(t;f)e df . (IV-1) 

, , 2 

The first important quantity is ( 6v(t) ), which, when compared 

2 

to |{v(t))| , is a measure of the intensity of the disturbance. The 
second-order moments <6v(t) 6v (t')) and (fiv(t) 6v(t^)) follow dii'ectly 
from the similar quantities for h(t,f) by applying Eq. (IV-1). We note, 
however, that from Figure IV-1, h(t,f) includes the effect of the trans- 
mitting and i-eceiving antennas, which are spatial filters, and the 
receiver, which is a temporal filter. 

To eliminate these detennlnistic factors we can compute the response 
of the transionospheric channel to a single monocliroraat ic plane wave for 
a time-invariant scattering layer. This is a moi’e fundamental quantity 
chan h(t,f), and we shall see in Section IV-D that h(t,f) can be readily 
derived from it. Thus, our problem x-educes to one of calculating the 
second-order complex moments of the random field that x'esults from pi'opo- 
gating a plane wave through a weakly scattering irx'egularity layer. 
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To formulate that problem, we first consider the equation for wave 
pi'opagation in an irregular medium — namely, 

e E = - • V log c) . (IV-2) 

The I’elative permittivity e is a random function parameterized by position 

and time, although in deriving Eq. (IV-2) we have assumed that the temporal 

variations of g are slow compared to the propagation time through the 

medium. We write e = (e) + 6si ^iid restrict ourselves to the class of 

2 

problems for which (e) s* 1, and <6e ) < < (e). (The angular brackets 

denote ensemble average.) It is then natural to consider a series solu- 
tion to Eq. (IV-2) with terms 0 in magnitude. 

The gradient term in Eq. (IV-2), which produces polaz'ization effects 
(Strohbelin and Clifford, 1967), is usually neglected on the grounds that 
significant changes in e occur only over large distances compared to a 
wavelength. For our applications, this is unnecessarily stringent. 
Following Balser (1957) and Tatarski (1971) , we have applied a series 
solution to Eq. (IV-2) directly. The result is a sequence of "cascading” 
inhomogeneous vector differential equations. 

If tlie Incident field (the zeroth order terra) is a plane wave 

—* 

E (r) = A exp f-ik • r j , it can be shown by repeated application of 
o o’- 

Green's theorem that the first-order terra E^Cr) is given by the integral 
expression 


E^(r) 




J (3)j6eCr ^)exp{-ik 


?'}G(r,?Od?' 


+ terms 0 



(IV-3) 


SO 



where 


GCr,?0 


exp{-il«|r - r^|} 



( IV-4) 


and a Is a milt vector in the direction of r. We see from Eq. ClV-3) 
r 

that contrihntions from the gradient term are ultimately negligible at 
points well removed from the scattering layer. 

We shall consider scattering only in a small cone of angles about 

a , so that the term in square brackets in Eq, (IV-3) is approximately 

'^^’e 

equal to A^, Moreover, foi' the fi’equencies of interest, 6e £2 - — ^ AN , 
where i* is the classical electi’on radius, and aN is the local departure 
of the electron density from its mean value. With these simplifications, 
Eq. (IV-3) becomes 


E fr) = - r PC3)r AN (rOE CrOGCr,r')dr' . (IV-5) 

1 e « u e o 

We note that Eq. ClV-5) is linear in E (r) . 

o 

The validity of this model depends on the accuracy with which 

— ♦ — t — < — . 

E^(r) represents the random component of E(r). Clearly, the approxi- 
mation E(r) =s E Cr) + E, (r) is valid only If the condition 
o 1 

<E^ * E *> < < ClV-6) 

'1 1 ' o 

is satisfield. It is coimiion expei’ience , however, that the first Born 
appi’oxijnation [Eq. (IV-5)] gives an accurate representation of the random 
component of E considerably beyond the few percent limit Imposed by 
Eq. CiV-G) . On the other hand, we know from theoi’y and experiment that 
multiple-scattering effects must ultimately be considered. Nonetheless, 
our results can still be applied by using Uscinski’s (1968) method. In 
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effect, the weak-scatter results are repeatedly integrated in a multi- 
layer model. 

It remains to assign a statistical structure to AN and compute 

e 

appropriate statistical descriptors for E^. We have used the transverse 

spectral decomposition of AN as was done recently by Lee and Harp (1069). 

e 

In Appendix A we have compactly summarized the main properties of the 
spectral decomposition of complex homogeneous fields. The power of this 
approach is demonstrated by an application to the free-space propagation 
of random fields. 

Before continuing, we note that the major difference between our 
own analysis and the extensive analyses of deWol- (1972), Tatarski (1971), 
Lee and Harp (1069), Strohbehn (1968), and others is that we have analyzed 
the quadrature components of the diffracted field I'ather than their ampli- 
tude and phase. To obtain amplitude and phase statistics we invoke a 
gaussian signal-statistics hypothesis. It h..^> -ens that amplitude statistics 
are simply derived, while phase statistics are quite difficult because 
the statistics are generally non-Rician (Section V-B) . 

By comparison, most researchers in computing phase statistics have 
either accepted the Rytov solution (cf. Tatarski, 1971; Strohbehn, 1968), 
or they have made the small-phase approximation of equating the signal 
phase with the phase -quadrature component (cf. Lee and Harp, 1969), When 
the small-phase approximation is valid, all results agree, including 
those derived under the gaussian hypothesis, 

B . The Transfer Function for a Weakly Scattering Irregularity Layer 

Following Lee and Harp (1969) and deWolf (1972) we assume that 
AN^(r) is statistically homogeneous (at least in the wide sense) in any 
transverse plane z = const. Hence, it admits a spectral decomposition 
of the form 
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AN^(r) = J Jexp{-iK: • p} dgCtciz) 


CIV-7) 


where r = p + a a, and a • p = 0. We shall refer to |Ck;z) as the 
z z 

— ♦ 

transverse Fourier spectrum of (r) , 

e 

“#• — ♦ ”4 •»+ 

We let k • r = • p + consider a scattering layer of 

thickness L centered at the origin of coordinates. Substituting Eq. £lV-7) 
into Eq, (IV-5) and changing the order of integration gives the result 


L/2 

E^(r) = - 3’^ JII d5(K;zOexp[-iK • 

“L/2 


„ j -»2 2)1/2 

*' exp^-ik Ap + Az f 


X 6Xp{lk^iz} JJ 

-TO ^ip + Az ^ 


exp[-i[('ic + *^p3}dAp dz' ClV-8) 


where Ap - “ p and = z' - z. The integration over Ap has been 

evaluated by deWolf (1972) as 


where 



exp{-ikg(K:;k^) | Az| } 
gClcjk^) 


CIV-9) 


gCicjk^) = 


[- 


Ck + 


— t 2 2 

k^) /k^ 
T 


1/2 


CIV-10) 


53 



By using this result, we can write Eq. (IV-4) in the form 


> •— i — + 1-^ 

E (r) = E (r) * (p;k) 
1 o 




where 




(IV-12) 


and 



L/2 exp{-ik gCK;k^) |A2| 3 

= J dg(K:;a')exp{.1lc Az] dz^ 

“ L/2 ^ s(k]\) 


(IV- 13) 


We see from Eqs. ClV-11) , CIV-12), and (IV-13) that E^(r) is derived from 

* — ► ^ » 

E (r) by scalar multiplication with a complex random function ;lj (oik) 
o 

that admits a transverse spectral decomposition [Eq. (IV-12)] of the 
form discussed in Appendix A . It is worth noting here that the only 
difference between the Born solution and the Rytov solution is that the 
diffracted field takes the form E(r) = E^(?)exp[ijj^(p;k)] for the Rytov 
solution. 

We shall characterize E (r)/E (r) [Eq. (IV-11)] by the two-frequency 

1 o 

correlation functions 


and 






(ty) 




(IV-14) 


(IV-15) 
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The wave vectors and differ only in magnitude (jk| = k = 2jtf/A). 

By using the orthogonal increments property of dg (Appendix A) , we can 
easily derive our main result: 

exp{“iK: * (IV-16) 

(2it) 

where 




Ca) -(p) 


L/2 

) " JJ oClc;z^,z")exp 

-L/2 




expj 

f 

pi 



1 

Az'l + k^^^g(K;k^^^]| 


]} 

X 



gl 

■pCo'; 



dz ‘ d’zJ' 


(IV-17) 


The upper signs in Eqs. (IV-16) and (IV-17) are used for the R, integral 


and the lower signs for the B integral. In Eq, (IV-17) §(k;;z',z'') is 

TZ 

the "transverse" spectral-density function for the irregularities. When 
z' = z" , §(K;z^,z'') is purely real, non-negative, and symmetric. 


Equation (IV-16) can be simplified by changing variables to 2x 
= z' + z''' and X| = z^ - z''. We assume that 


§(1 c;z‘',z") = M-Cx) iCks'II) 


(IV-IS) 
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where 


|aCX) 21 


nCO) 

0 


|x| ^ L 
lx[ > ^ 


CIV-19) 


The result is 


= ^( 0 ) 


exp 


■| 2 izH (ll; 


k'“' 


fCci') ii(p) 


g(K;k^P^) 


J S (K:;Tl)exp |-i'nH'*‘^lC;k^‘^\k*'^^)|’slnc ^CL “ ^ 


CL - T|)d'i] 


n “* T */“♦ ^ 

+ J 9CKtT|)exp j -iTlH ^tc;k ,k 


Cp) 


sine 


;a + H)H(K:SW.k‘8> j 


(L + 'i])d'I] 
(lV-20) 


where 


H 


= / 2 , H* = [g'"’ ± g'P^ /2 


and 


(IV- 21 ), (IV- 22 ) 


Cev) Cev) 

, Cp) 

G ^ s= k ^ cos 0 


(cv) 


1 




(IV-23) 
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In the coordinate system of Figure IV-2 


s' 


sec 0 g 



2 2 

K + K 
X y 


(k cos 0 ) 



k cos 0 


( IV-24) 





FIGURE lV-2 GEOMETRY FOR LAYERED-MEDIUM SCATTERING CALCULATIONS 


A further simplification can be I'ealized when fCK:;T]) Si 0 for 

T1 > Tl < < L. We can then neglect the T|“dependence in the sinc-function 
“max 

terras in Eq, (IV-20) and extend the limits of integration to infinity. 

The integral in Eq. (IV-20) then simplifies to 
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(IV-25) 



-a* 1 yCo') 7*CB)\ 
f I k;1c ,k 1 

\ / 

\ / 


where 

00 

§Ck, IC^) = J |(KJT[) exp[i'r]lC^}dT| . (IV-26) 

“XO 


We note that in general the entire z dependence is contained in the 
first term of Eq. (IV-20) * 


To obtain some feeling for the physical meaning of these results 

consider the single -frequency case in which = f . Then, 

H = 0 for the R,. integral, and it is independent of z. This result is 
Vz 

expected because of the free-space results of Appendix A. With a suitable 

normalization, the correlation functions for the in-phase and phase- 

quadrature components are — [ l + RefB, 1], respectively. In this formu- 

2 ifz 

lation there is never a pure phase perturbation because propagation, 
which takes place continuously, acts to redistribute the phase perturbation 
as they occur [see Appendix A], 


In the next subsection we shall specialize these results by intro- 
ducing a specific spectral-density function appropriate to the ionosphere. 
We shall then present some simplifying approximations that are commonly 
employed and discuss their validity. 


C . The Special Case of a Power-Law Spectral-Density Function 

Satellite in-sit u electron-density measurements reported by Dyson, 
McClure, and Hanson (1973) show that a power-law spectral-density function 
is a reasonably accurate analytic form for §. This is supported by the 
analysis of radio-star scintillation data by Cronyn (1970) and 
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Rufenach (1972), We have chosen (after Tatarski, 1971) the analytic form 


5(S ) = Stt 


3 / r(Y + 1/2) \ 

' ./^rcv “ f 


./tTFCy “ 


[1 + ((yS) ] 


Y + 1/2 


(IV-27) 


where 


= f B'(Y)tC^ + tC^ 1 + / B(Y) - 2k K C(Y) 
\ X y / 2 X z 


(IV-28) 


2 2 
B^(Y) = a cos Y -3- sin Y 


(IV-29) 


and 


2 

C(Y) = (1 - a ) cos Y sin Y 


2 2 2 

B-(Y) = cos Y + a sin Y 


(IV-30) 


(IV-31) 


The angle Y is the angle to the principal irregularity axis. In 
the ionosphere it is the geomagnetic dip angle. The parameter a is the 
axial ratio (a & 1), The parameters o' and Y irregularity scale 

and spectral index respectively. The spectral-density function [Eq. (IV-27)] 
is normalized so that 


^ dK dK dK 

! (3) J i(S-) 17 17 " 1 • av-32) 

The anisotropy was introduced, following Sudden (1965) , by performing a 
coordinate transformation of the isotropic spectral-density function 
(see Appendix D) . 
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By substituting Eq. (IV-27) into Eqs* (IV-25) and (IV-16) and 

changing variables to k ^ and k £v we obtain the result 

X y 


R 


B 


= ± cr^ sec 0 


[s/ ( rOL±^ 
L \/«r(Y“i) 


)•■] 




exp [-K . Ap] 


Ck:)<p 


dK 


( 2 jt) 


(IV-33) 


where 


and 


= 1^(0) r^ sec 0) <AN^> h 

i B ' e' 


5 = exp 


2iz(k cos 0) 


CIV-34, 


( 1 + l\ 
— 2 — ) 


Sine { ^ k cos 0 h) 


^ r 2 ,V + V2 , 

[l + S ] 


(IV-35) 


In the new coordinates, 


2 / 2 2 \ 2 

S = + Kyj + [H*p] BCV) - 2K^H*pC(Y) ClV-36) 
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and 


(a) 


1 “ 


2 2 

K. + K 
X y 

2 _ 2 
P (1 + Ak/k) 


2 tan 0 


cos cp + Ky sin cp) 


p(l + Ak/k) 


1/2 


(IV'^37) 


where p = o' k cos. 0 , The parameters are summarized in Table IV-1. 

2 2 

The normalization factor h is chosen so that E n /2Z is the 

o T o 

scattered power per unit area crossing a plane perpendicular to z [cf* 
Rino and Fremouw (1973)]. The linear K ,K -term in Eq. ClV-35) was 

^ y 

introduced to reference the calculations to the ”line-of— sight” intercept 
in the z plane (see Figure IV-2) , which is the natural reference point. 


The parameter p is the ratio of the scale parameter o’ to the pro- 
jected wavelength X sec 9 , If it is sufficiently large, we can apply 
the so-called Fresnel approximation 


1 


(a) 

/Cp) 

- g' ^ ^ 


2 2 2 
k: + iC + tan eClC cos CP + 

X y X ^ 


2p^(l + Ak/k) 


2 

K sin cp) 

y ^ 


tan e Ck; cos (p + fC sin qj) 

+ . (IV-38) 

p(l + Ak/k) 


By using Eq. (IV-38) in Eq. (lV-23) and the equivalence Ak/k = Af/f, we 
obtain the approximate forms 
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Table IV“1 


PARAMETERS FOR COMPLEX-CORRELATIOK-FUNCTION FORMULAS 


Parameters 

Symbol 

Comments 

Wavelength Cs) 

, , (o') , (p) 
k;k ^ 

k^°"^ = k - Ah, k + Ah 

^Note that Y " 

RMS electron 
density 

. 2 , 


Distance 

z 

Distance from ionosphere 
penetration point to re- 
ceiver, or corrected dis- 
tance for satellites 

Axial ratio 

a 

Anisotropy parameter (a ^ 1) 

Magnetic-field angle 

Y 

Geomagnetic dip angle (a- 
ionosphere penetration 
point) 

Scale size to wave- 
length parameter 

P 

p i O' h cos 0 , where a is 
power-law spectral -density 
function scale factor 

Spectral index 

Y 

(Y > 0) 

Incidence angles 

9- 

CD 

Theta, 0, is the zenith angle 
at the ionosphere penetra- 
tion point. Phi, qj, is the 
azimuth measured from the 
magnetic-meridian plane. 

Layer thickness 

L 

Centered on origin of 
coordinates 


G2 







H =- 


r 




P 9 O 

K + k: + tan^e 

X y 


Ck 

X 


2 

cos cp + ICy sin cp) 


2 


2 

P 



Af/f 

1-CAf/f)^ 

1 


for R 

'I' 


(IV-39) 


4 tan 0 



cos CD 4 1C 

^ y 


p 


sin cp 


for B 


For H4', the opposite values are used for R and B [cf. Eqs. (IV-21) and 
(IV-22)], 


When these approximations are valid, we obtain the "large p" form 
of Eq. CIV-35), 


r t- 


4 2 : exp 

z 


A 


iZ 


Af/f 


1 - (Af/f)' 
1 

1 - CAf/f)‘ 




2 2 2 2 
K 4 K 4 tan 0(tc cos cp + K sin m) 

X y X ^ y ^ 


X sim- 


- Ck cos 0 ) 




[l 4^B^(Y)K^4ic^)4 (H*p)^BCY) -2 h;^(H*p)CC'F)] 


where 


Y + i/2 
(IV-40) 


Z = 


Az sec B 

„ 2 
2jr O' 


(IV-41) 
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In Eq. (IV-40) H*p is negligible [see Eq. (IV-39)3 for the Integral, 

and it is equal to tan 0 Ck cos m + ic sin (n) for the R integral. The 

X ^ y ^ ^ 

effect of the H*p terms is to distort the angular spectral-density 

1 

function which, from Eq. (IV-35) , is given as 7 -^ . 

9,r Y + 1/2 

gll + s"=] ' 

r 2 -] y 1/2 

In Figure IV-3 we have plotted log [1 + S J' versus log 

2 

(K /p) j with K = Slid S given by Eq. (IV-36) with a = 5 and Y = 45°. 

X ^ ’ y 

We set y = 1.5, which is close to the value inferred from satellite 
measurements (Section IV-D) . The straight lines have a slope of 4, which 
is the assnnptote when the p-terms are neglected. We note that for the 

smaller p values, the departure from the asymptote is significant when 

, 2 
(t/pjis large. Numerical calculations show that 0 ^ is underestimated by 

approximately 50% when p = 1, but it is within a few percent of its 

correct value when p = 10. Evidently Eq. (lV-40) is valid when p is 

greater than 10. Hence, in the remainder of our calculations the 

simplified form of Eq. (IV-35) will be used. 


CN 

Ui 

4* 


O 

cT 

o 

-1 

+ 


LOG^q (Ky/p) 



FIGURE IV-3 SPECTRAL-DENSITY-FUNCTION CROSS-SECTION PLOTS 
SHOWING THE EFFECT OF H*p TERMS IN EQ. (40> 
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D. 


Effects of Time-Varying Media 


In almost all applications we measure time series that arise because 
the medium or the source or both are moving relative to the receiver. 
Moreover, the medium undergoes internal rearrangement, although we can 
usually neglect this effect. In Appendix B we have derived general 
formulas for converting two-dimensional spatial results to one-dimensional 
temporal spectra. 

As an example, consider an in-situ measurement with a moving probe 

such that v =0, The power spectrum derived from Eq. (IV-27) is then 
z 

proportional to 


1 + 


2 2 
4itQ< ^ 


U 


Y - 1/2 


1 - 


X 


2 2 

u + B^(Y)u 
X y; 


(IV-42) 


The effective velocity is 



2 


u 


X 


X y 



2 “ 

which is approximately equal to u when B'(Y)u > > u • However, since 

y* y 

B^(¥) > > 1 for highly elongated irregularities when 7 os not too large, 

the condition B'(w) 1 can be satisfied even if u and u are comparable 

y X 

in magnitude. 


In a transmission experiment the potential observable is the complex 
voltage measured effectively at the antenna terminals. Since the angular 
width of the field spectrum is typically much narrower than the antenna 
beamwidth, however, the complex correlation functions are proportional 
to the Ap = 0 terms in Eq. CIV-33) or the 6 x = 6 y = 0 terms in Eq. (B-4) 
of Appendix B. Spaced-receiver measurements then give the dependence 


directly. 
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Unfortunately, the corresponding complex poiwer (temporal) spectrum 

Integrals for Eqs. (IV-35) or (IV-40) cannot be evaluated analytically. 

However, for normal incidence one can show that § (^) is proportional to 

z 


A /4? 



X correction factor. (IV-43) 

In Eqs. (IV-42) and (IV-43) we have used 'I?*’' for temporal frequency to 
avoid confusion with the RF frequency f. The correction factor will 
exhibit a weak temporal-frequency (f) dependence. 

We first note that the effective J'-index in Eq. (IV-43) differs by 

r — (2V — 1) . — 2 yT 

unity from that inferred from Eq. (IV-42) \ j vb‘. J J. This 

has been pointed out by Cronyn (1970). Dyson et al. (1973) foxrad an 
-1.9 

j frequency dependence from their in-situ measurements. Hence, 
from Eq. (IV-42), y = 1*4S, which is one-half the index we would 
observe by Fourier-analyzing the temporal fluctuations of the diffracted 
signal. One must keep in mind, however, that there is no simple relation 
such as Eq. (IV-42) for interpreting complex temporal-frequency spectra 
that include the Z-dependent terras. 
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On the other hand, the first-order moments are vmaffected by the 
"scan" because they are proportional to integrals over all temporal 
frequencies. It is easy to see from Eq. (B-4) that the results are 
Identical to those obtained by integrating over all spatial frequencies. 
Hence , the formulas 


sec e I [1 t Re[B^(0;f)]l 


(IV-44) 


and 


2 2 
Oj sec 0 '* 


C 

xy 


1 2 
2 


sec 0 JmfB (0;f)} 
'll 


(IV-45) 


(IV-46) 


■which were used in Section III correctly predict the variance of the 
in-phase and phase quadrature components of the signal and their 
covariance . 


E. Method of Humeri cal Computation 

The results presented in Section IV-C give the spatial and radio- 
frequency or wavelength dependences of the diffracted field from a 
weakly irregularity layer. The spatial dependence is derived from a two- 
dimensional Fourier transform [Eq. (IV-33)]. However, in the remainder 
of this report we shall consider only the Ap = 0 term. In Section IV-C 
we showed that these integrals are generally adequate for single-receiver 
measurements. 

Hence, the results that we shall be interested in demand an inte- 

(o) 

gration of a [Ea. (IV-35)] over all spatial frequencies for which g' 

(a) ^ 

and g^ ^ [Eq. (IV-37)] is real. Note that, in general, four integrations 
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must be performed to obtain RefR^^j}, jftii[R^], Re[B^|j}, and . Moreover, 

whenever the Z term is finite (that is, when H 0), the integrand 
oscillates. Indeed, the integrals are of the most difficult type to 
evaluate . 

The large p approximations are of some help, but they do not elimi- 
nate the oscillation problem. Integrating over the polar coordinates 
5= k: cos s and = k sin o, however, does simplify the integration 
considerably. The formulas that we have programmed for machine inte- 
gration can be summarized as 


where 




KF(K:)dlC 



(IV-47) 




pF 1 

' 

FCk) = r exp 


R 

2 2 2 
k: [ 1 + tan 0 cos (® - g^)] 



F 

B_ 



X sine ) 


pF 1 
R 

P 

2 2 2 / 1 
1C [l + tan e cos ($-gi)] +K\~ 

- Jtan 0 cos (* - |2f) 

Pl?l] 

2 


j. 

B 

\ ' 


L ^ J 


2 / 2 2 \ 2 2 
1 + K ^B'(Y)cos f + sin s) + [H*p] B('i')-[H*p] K:cos*C(y) 


-Cy + 1/2) 


d6 

( IV-48) 


£ 


XL sec 9 
2 

23TQ' 


CIV-49) 
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(IV-50) 


F = 

^ 1 - 


F = ^ av~51) 

1 - (Af/f) 

tan 0 cos(q - pf) for R 

(IV-52) 

0 for B 

No further simplification is justified unless the results are 
restricted to normal incidence. In that special case both the 2 and 
£“dependent terms in Eq. (IV-4S) are §-independent . Hence, the §- 
integration is simplified considerably. The argument is a smooth, 
slowly varying function of ®, and almost any numerical integration scheme 
will work. To evaluate the Z-dependent integral it was necessary to use 
an adaptive Simpson rule (see for example, Lyness, 1970). Basically, 
the adaptive methods vary the mesh size according to how fast the inte- 
grand is changing. 

Finally, if the £-dependent term is ignored, and the restriction 
to normal incidence is retained, it is possible to obtain a series 
solution involving Legendre functions and ellptic integrals that can 
be evaluated with sufficient accuracy by approximate formulas. The 
results, which are easily programmed for machine computations, ax’e 
summarized in Appendix C. 

These results show that the most troublesome factors result from 
oblique incidence, which, unfortunately, cannot be ignored. Rino and 
Fremouw (1973) presented a general formula that allowed oblique 
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incidence. However, in that foiTuula the k;^ term of the spectral density 
function was set equal to aero. We see from our present results that 
this is valid only for normal incidence, when the [H*p] terms can be 
neglected. 

However, the error affects only the integral [cf. Eq , (IV-52)], 

so the formula for B given by Rino and Premouw (1973) is correct. The 

[H*p] terms affect through m. in Eq. CIV~34) . Hence, the incidence- 

2 

angle dependence of cr^ is in general somewhat more complicated than a 

sec 0 variation. It could be important, for example, if one were to 

2 

predict or evaluate the variation for an orbiting satellite. 

In the next subsection we present representative results derived 
by numerically integrating Eq. (IV-47). Some additional calculations 
will be presented in Section V. 


Discussion and Examples 


To characterize the statistics of a plane-wave field diffracted by 
a weakly scattering irregularity layer, we have computed the complex 

covariance functions ^“'^f j and j as 

defined by Eqs. (IV-14) and (IV-IS). The Ap dependence is given by a 
two-dimensional Fourier transform [Eq, (IV-16)], and we have considered 
only the Ap = 0 correlations. 

The general results depend explicitly on each pai'ameter in Table IV-1, 


However, the approximate results [Eqs. (IV-40) , ClV-41)] show that the 

2 Xz sec 9 

TT 

2 


/ 2 Xz sec 0 \ 

ratio of the Fi-esnel-zone area I I to the transverse scale 


parameter q» , which we have denoted by Z , is fundamental. Its value 
determines the degree to which propagation has altered the structure of 
the field from what it was in the vicinity of the scattering layer. 

Indeed, a general chaz'acterization of the statistics can be made 
in terms of Z. Consider that, independent of all other parameters. 
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Re[B^3 negative and near or equal to its maximum absolute value of 
unity when Z - a. Moreover, = 0 when Z - 0. Now, 


as Z increases from zero to infinity, RefB } tends to zero, while 

'I' 

epra[B 3 first becomes more negative and then tends to zero. 


The behavior of Re{R^3 depends on AS (we still assume 


that 


Ap = 0). When A^ = 0» is Independent of Z and equal to 

unity. When A^ is finite, Re{R^} decreases from unity to zero as Z 
varies from zero to infinity. Similarly, j?m[R^} is zero when A^ = 0 
and for finite Af it varies from zero to a positive maximum and then 
decays to zero again as Z varies from zero to infinity. 

To apply these results, we let the diffracted field at some 
frequency f be represented by E and E' at some other frequency f 
It then follows from Eq, (IV“33D that 


1 2 


R^= <XX'> - <X> (XO = 2^T } + Re{B ]] (IV-53) 


R 


1 2 


yy 


= <YY'> - <Y> (Y^) = 2 f^T ClV-54) 


R = <XY'> 
xy ' ' 


“ <X) <Y'> = 


1 2 

r c? eWR,} -^tn^B }] 

* 1 y 


(IV-55) 


and 


R = (X^Y> 
yx 


- (X') <Y> = “ ^ sec eDmfBj +^m{B, j] (IV-56) 


2 "T 


where X and Y denote the real and imaginary parts of the field. These 
formulas are simply generalizations of Eqs. (IV-44) , (IV-45), and (IV-46), 
Indeed, when f = f', = a®, = dy, and = C^. 
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Now, from the limiting behavior we have just described, it follows 

that when E is sufficiently small, R > > R , and both R and R 

yy XX g xy yx 

are near zero. When f = f', this condition becomes n > > cj®: that is, 

y X 

the scattered power is nearly in phase quadrature with the undeviated 

component. In Rino and Fremouw (1973) this condition characterized the 

"near zone." When Z is sufficiently large, R and R are nearly 

yy XX 

equal, although when ^f is finite the limiting value becomes smaller 

2 2 12 

and smaller as Z increases. When = 0, q sec Q. This 

X y 2 T 

behavior characterized the "far zone" in Rino and Fremouw (1873) , 

All parameters in our channel model (Table IV-1) that can be varied 
independently of Z will only affect the rate at which we move from near-’ 
zone scattering to far-zone scattering as Z is increased from zero, 
to infinity. We have illustrated this below by using the equation 
presented in Section IV-E. 

1 . Normal- Incidence Computations with Af - 0 

To evaluate Eq. (IV-48), we must specify the axial ratio a, 
the dip angle Y, a layer thicloiess parameter £, the spectral index y, 
and the incidence angles 0, cp as well as [Eqs. (IV-47) through 
(IV-52)], Wlien 0 is finite, it is also necessary to specify the ratio 
L/2q> [Eq, (IV-48)]. However, to start with we shall consider only the 
normal-incidence single-frequency cori’elations for which Af = 0 = cp = 0- 

In Figure IV-4 we have illustrated the effect of the "layer- 

smearing" £-dependent factor in Eq. (IV-48) . The dashed line is for 

zero layer thickness, and the solid line is for conditions we expect 

to be typical of ionospheric scattering. The effect of finite layer 

thiclcness is confined to the near zone. It acts primarily to make 
2 

cr , finite when Z = 0. 
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FIGURE IV-4 

In Figure IV-5 we sliow the effect of varying the axial ratio a . 

We see that as a increases, the convergence to far-zone conditions 

becomes slower. We expect this, since increasing a effectively compresses 

the spectrum toward the lower spatial frequencies. Hence, Z must be 

comparatively larger before propagation can have an impact. Evidently 

we could produce the same effect by increasing the spectral index y. 

2 

This is verified in Figure IV-G, where we have plotted 0 as a function 

3C 

of Z for different values of the spectral index y Varying the Incidence 
angles (9, qj) has a similar effect, although the details are more compli- 
cated. Because of the increased complexity of the computation, however, 
we have not presented numerical results for Oblique incidence, 

2. Coherence Bandwidth 

If two monochromatic signals at frequencies f and f ' are 
transmitted through the ionosphere and the received complex signals 
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FIGURE lV-5 FAMILY OF vs. 2 CURVES SHOWING EFFECT OF ANISOTROPY 



LA-2273-23 


FIGURE IV-B FAMILY OF vs. Z CURVES SHOWING EFFECT OF 
POWER-LAW SPECTRAL INDEX 
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are cross-correlated, one measure of coherence bandwidth is given by 
I (E [ . We do not choose to subtract the mean signal values, since 

one normally processes the complete signal. 

From Eq. (IV-33) it follows that 


(E E'*) 


(E> <E^*> + sec fi (Af,f) 


(IV-57) 


However, if Af is not too large, from Eq. (IV-34) it follows that 


(E E^*> 


25 1 + cr„ 


sec GfR (Af,f) 

It 



(IV-58) 


Hence, we see that coherence bandwidth cannot be specified independent 
of the scattering intensity. 

Nonetheless, the At frequency dependence is contained in the 

covariance function R, . In Figure IV-7 we have plotted 

^ z 

function of A^/f for two different values of Z. As we should expect 
from the discussion at the beginning of this subsection, the results 
depend critically on the value of Z. However, since Z itself is fre- 
quency-dependent, the curves in Figure IV-7 are not universal. As the 
center frequency, f, is increased, Z decreases if all other parameters 
are held constant, 

2 

The implication is that for a given value of the coherence 

bandwidth will be larger at a higher frequency than at a lower frequency. 
In Section V we shall calculate the two-frequency correlation function 
for intensity for comparison. 




as a 
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V APPLICATIONS 


In Section IV we derived general formulas for the complex correla- 
tions that characterize a plane-wave field diffracted by a weakly scatter- 
ing irregularity layer. We first calculated the time-invariant correla- 
tions <E (^;f)E *(p;f')> and <E (^;f)E (p^;f')). We then considered the 
z z z z 

generalizations that allow for time variations (see Appendix B) . 

Hence j the most general observable is a time function with position 

^ ^ 

and frequency dependences— say , V (p;f) = X Cp;f) + iY (p;f). For 

t t t 

notational simplicity we have suppressed the distance parameter z. The 
position and frequency dependences can, in principle, be verified by 
spaced- receiver measurements with multifrequency transmissions. However, 
in most experiments only amplitude or intensity is available for each 
independent position-frequency measurement. For this reason, intensity 
statistics are most important for channel diagnostics. 

In this section we shall apply the results of Section IV together 
with the gaussian signal-statistics hypothesis to derive the second-order 
moments of intensity. We shall also discuss the method of making a similar 
computation for phase statistics. 

A. Intensity Statistics for Gaussian Fields 

The observable of interest here is the intensity of a complex gaussian 
random field. Note that in the theory of stochastic processes, a random 
field is simply a function parameterized by one or more variables. The 
height of an ocean surface is an example. In our case we have a complex 
function of position, frequency, and time. The power of the gaussian 
hypothesis lies with the fact that the joint statistics of X and Y 
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depend only on the complex correlation functions calculated in Section IV 
or their generalizations derived in Appendix B. 

In what follows we assume that the observed field is gaussian and 

2 2 2 

normalized to unity average power. Recall that o “ then be 

interpreted as the fraction of the incident power that is scattered. 

We let I & W* . A primed dependent variable will imply that it is evalua- 
ted with the independent variables p , f , and t , In Appendix E we show 
that 

<II0 “ <I>^ ^ 2a^Cl - 0^)[Re{R} + 1 b|cos 2(^ - p{)] 

+ + |b|^] (V~1) 

where R and B are given by Eq. (IV-33) or Eq. (B-4). We have also used 
the definitions tan 2Q ss J3ra[B}/Re{B] and tan = (Y)/(X) . We note that 
for interpreting most experiments, ^ can be set equal to zero. That is, 
the coherent component is the phase reference. 

When eJ is equal to zero, the term in square brackets reduces to the 

correlation of the real part of V, R [cf, Eq. (1V~53>]. Hence, when 

2 2 
Q « 1, the quadratic terms in q can be neglected, and Eq. (V-1) 

reduces to 

<IlO " <I>^ = 2 ct^R . (V-2) 

This result can be derived directly by equating the amplitude perturbation 
with the in-phase signal component. 

The Fourier transform of Eq. (V-1) or Eq. (V-2) gives (in expectation) 
the power spectrum of the intensity fluctuation if I = I', or the co- 
spectrum if spaced receiver and/or two-frequency correlations are employed. 
The approximate form given in Eq. (V-2) has been used extensively for data 


78 



intei’pretation [see for example, Cronyn (1970) and Rufenach (1972)], 
although even then the interpretation is not simple [see Section IV-D]. 

We also note that when I = l', Eq. (V-1) reduces to the formula 
for the scintillation index- Indeed, for that special case, R s 1, 
and Eq. (V-l) becomes identical to the formula given in Rino and Fremouw 
(1973) that was used in Section III [cf. Eq. (III~11)]. In the simplest 
case, Eq, (V-1) gives the temporal structure of the signal fades whose 
statistics we analyzed in Section III. The formula Is, however, more 
general since it gives the spatial and frequency dependencies as well. 

To summarize, in Section III we discussed the method of computing the 
first-order amplitude statistics of a complex field. In Section IV we 
computed the complex second moments of the field. Finally, in this 
section we have made use of the fundamental propei'ty of gauss ian fields — 
namely, that the complete statistical structure of the field can be 
determined from the complex covariances R and B. 

Before presenting some examples derived from Eq , (V-1), we emphasize 
that in general, R and B depend onAp = p - p\ f = [f + f']/2, Af = f - f', 
and 6t = t - t'. In general, there is a more con^ilicated dependence on 
6t which we have ignored by assuming "frozen fields" [Appendix B] , When 
Ap = Af = 0, the Fourier transform of Eq. (V-1) gives the power spectrum 
of the intensity fluctuations. The width of this power spectrum is a use- 
ful measure of the fading rate, which is fundamental for designing receivers. 

Unfortunately, wliile the statistics of any measurable quantity can 
in principle be computed from R and B, other results are not as simple 
as Eq. (V-1). Phase statistics will be discussed in Section V-B. 

1- The Wavelength (Frequency) Dependence of S^ 

— — ^ 

Since the S^ scintillation index is the most commonly used 
4 

measure of scintillation activity, its frequency dependence is easily 
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measured. We first note that a « X . Hence, in the S scintillation- 

2 2 2 ^ 2 22 
index formula it is convenient to replace a by \ where a* 4 cf /X • 

,2 

We can think of a as the coherence ratio for X = 1 m. Recall that for 

2 

unity incident intensity, a can be interpreted as the fraction of the 
incident energy (or power) that is randomized. 

Similarly, from Eqs. (IV-41) and (IV-49) , we see that Z and £ 
are proportional to X. Hence, they contain the entire wavelength 
dependence of B , and therefore, the remaining wavelength dependence of 

'll 

S [R = 1 when Af = Ol. We define Z* = Z/x and £ = £A in a manner 
4 i[r 

analogous to the definition of 


Ill Figure V-1 we have plotted the logarithm of against 

,2 

the logarithm of \ for different values of cr' for a power-law spectral- 

density function with y = 1.5. The curves show a decreasing slope with 

increasing S. . Moreover, the results show the rate of decrease of the 

2 

slope to be largest for the largest values of a . Ultimately, the slope 
becomes negative. However, the theory has not been validated for strong 
scattering conditions . 


Aarons, Allen, and Elkins (1967) define a spectral index for 


scintillation as 

iog(s^/sM 

4 4 

^ log(XA^) 


(V-3) 


Actually, their definition is in terras of frequency rather than wave- 
length, which changes the sign of T|. That is, Eq. (V-3) implies a rela- 

Ti -n 

tion of the form S « X or S « f depending on whether frequency or 
wavelength is used. 

Aarons, Allen, and Elkins (1967) report the results of radio- 
star scintillation observations at four frequencies spaced from HF to UHF. 
They I'eport that Tl decreases from ~2 for small index values through zero 


SO 




FIGURE V-1 WAVELENGTH DEPENDENCE OF S 4 SCINTILLATION INDEX 

for large values just as the curves in Figure V-1 predict. It should 

be pointed out, however, that extrapolating the results to L and S bands 

does not account for the indices that have been observed at those 

4 

frequencies. Evidently, the spectral shape deviates from the power-law 
model in the small-structure region. 
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2 , 


Two-Frequency Correlatioa Function for Intensity 


As a final example, we have computed from Eq, (V-1) the two- 
frequency measurement correlation function for intensity. The results 
are shown in Figures V-2(a) and (b) . The curves are normalized to unity 
at Af = 0. 

In Figure V-2(a) we show curves corresponding to the pt.'-ameters 

2 

used in Figure IV-7. These curves were computed with u =0.25. We 

see the same general property — namely, that the signals decorrelate 

faster for larger Z values. Moreover, E decorrelates faster than i. 

2 

In Figure V-2(b) we have repeated the computation with cr = 0.5 to show 
that the intensity decorrelate faster as the level of activity increases. 

B. The Method of Computing Phase Statistics 

To conclude this section, we discuss the method of computing phase 
statistics. Consider a complex gaussian random process evaluat'd at two 
"points" giving X = (X, Y, X' , Y^) as a data vector. The covariance 
matrix for g is 

A = “ </>) (x - <^>|> = 


where T denotes transpose. The joint probability-density function for 
X, Y, , and y' can be written as 

P^(X,Y.X',Y0 = “ <X^>)a"^(x - <X>)j- . (V-5) 
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FIGURE V-2 TWO-FREQUENCY CORRELATION FUNCTION FOR INTENSITY 
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These results are well known; they are discussed in any basic 
statistics textbook. By a simple change of variables , it is easy to show 
that 


PCA,?J,A' = A A'P CA cos 9 S, A sin izf, A'cos lZ^',A^sin gf' ) . <V-6) 

The joint probability-density function for ^ and eJ' is then given formally 
as 

CO CD 

PCjrf.pt') = J J P(A,KJ,A',pfOdAdA' . (V-7) 

o o 

Finally, from Eq. (V-7) it follows that 

=//!!{ pf'P(pJ,£f')dgfd)2J' . (V-8) 

For the special case of Hician statistics, it is possible to obtain 
an analytic form for [See, for example, Davenport and Root, 1958, 

pp. 161-165,] However, we have shown that one cannot assume that the 
quadrature signal components are of equal intensity and uncorrelated. 

At the present time iL is not known whether similar results can be 
obtained for ger.c:ral gaussian statistics. The formal results are pre- 
sented here mainly to illustrate further the technique of combining the 
gaussian-statistics hypothesis with the results of Section V to calculate 
various measurable quantities. 
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VI SUMMARY, DISCUSSION, AND RECOMMENDATIONS 


Through the research conducted under this contract we have (1) 
evaluated a gaussian signal-statistics hypothesis for ionospheric 
scintillation, (2) reformulated the diffraction theory in terms of 
quantities that can be employed directly in the analysis of linear 
systems, and <3) presented examples that illustrate the characteristics 
of the trans ionospheric channel. 

A. Gaussian Signal Statistics 

The main thrust of our work has been to characterize the structure 
of the complex random field that results a a wave propagates through a 
randomly irregular medium. The most natural and useful a-priori assumption 
is that the quadrature components of the complex field have jointly 
gaussian statistics. Historically, gaussian signal statistics were indeed 
the first to be postulated (Section II). 

Since only amplitude or intensity data are available in most experi- 
ments, however, it is usually not possible to test the gaussian hypothesis 
directly. In Section III we presented a formula for the second central 
moment of intensity and showed that it depend on quantities that could 
not be unambiguously determined from intensity data alone. This problem 
does not arise for the special cases of Rayleigh and Rice statistics. 

However, researchers in radio astronomy and optics have consistently 
observed that Rayleigh and Rice statistics do not correspond to measured 
intensity histograms. This fact alone is not sufficient for rejecting 
gaussian statistics, but there is a considerable amount of theoretical 
and experimental evidence supporting the competitive log-normal hypothesis. 
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According to the latter hypothesis, the logarithm of the complex field 
has jointly gaussian statistics for its quadrature components rather than 
the complex field itself. 


Because of our earlier work [summarized in Rino and Fremomw (1973)] 
and the work of Armstrong, Coles, and Riekett (1972), we believed that 
the differences between scintillation-data histograms for gaussian and 
log-normal statistics would be subtle and difficult to ascertain. Thus, 
one of the main tasks was to carefully evaluate intensity histograms. 

The results presented in Section IV-D show that with the one exception 
out of 10 cases, the best-fit gaussian PDF achieved a lower value for its 
chi-square-fit parameter than did the corresponding log-normal PDF. The 
exceptional data set can possibly be discounted because of its incon- 
sistency with simultaneous data taken with a separate receiver (Section 
III-D-2). In addition, for a given scintillation intensity, as measured 
by S^, for example, interplanetary and ionospheric scintillation data 
show a striking degree of similarity in their intensity histograms. 

Unfortunately, the data base is not as large as we had originally 
planned, as explained in Sections III-D-1 and III-D-2, Nonetheless, 
the available data do support the following conclusions; 

• Sc inti Hat ion- intensity probability-density functions 

derived under the assumption of joint gaussian statistics 
for the quadrature-field components achieve a better data 
fit than the corresponding log-normal probability-density 
function. 

■ The signal structure is characterized by the fact that , 
independent of the scintillation level (e.g. , S^) , more 
than SO percent of the random signal component is in phase 
quadrature with the steady or coherent component. An 
equivalent statement is that the equiprohability ellipses 
characterizing the random component of the signal have 
their semi-major axes more than ten times as large as their 
semi-minor axes, and their orientation is within 10° of per- 
pendicularity to the coherent signal component. 
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■ The log-normal probaTaility-density function tends to 

achieve a smaller rise time than the measured histograms and it 
underestimates the peak. In the tail region, the differences 
are not significant. 

• The signal structure that we have described is independent 
of incidence angles and the detailed shape of the spectral- 
density function that characterizes the irregularity medium, 
provided that it is smooth. However, structure comparable 
to the Fresnel-zone area must be present. 

These results must be qualified by the fact that no reliable data 
were processed with indices greater than 0,55. The limitations imposed 
by the implicit first Born assumption will be discussed in ’detail in 
Section VI-D below. 

B . The Hakagami Distribution 

Because the Nakagami distribution has been used extensively in 
ionospheric scintillation studies, its applicability needs to be discussed. 
It is usuallv described as an approximation to a general gaussian pro- 
bability-density function for amplitude or intensity if the appropriate 

variable change is made. It is convenient in that it depends on the 

-2 

single parameter m = S . 

We have pointed out, however, that the exact relationship between 
the Nakagami distribution and the more general gaussian family of dis- 
tributions is not established. For example, one cannot specify the 
joint statistics for the quadrature signal components that will give 
rise to a Nakagami amplitude distribution. Because of this uncertainty, 
the Nakagami distribution should be viewed simply as a convenient dis- 
tribution that can be easily applied to scintillation data. 

As such, its applicability depends mainly on how well it works for 
the intended application. We showed in Section III-E that, like the 
log-normal distribution, it tends to be broader than the correct gaussian 
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distribution. We have made direct comparisons for fade-margin deter- 
mination (Section III-E) to show that the ITakagami distribution gives a 
conservative bound for the fade margin. 

For systems planning, a conservative margin is probably desirable. 
Moreover, the differences are small enough (less than 1 dB at the 90 
percent margin for S s 0.5) to be of no practical consequence. We 
suspect that a similar conclusion would hold for scintillation-index 

conversion, at least for the smaller values of S,. 

4 

C . Complex Second Moments for Ionospheric Transfer 

Function Characterization 

Basic to any channel-modeling effort is a means of calculating the 
diffracted field for an elementary incident wave. We first considered 
a stationary time-invariant medium and then generalized the results. 

We have chosen to compute the complex second moments of the diffracted 
field, which are most convenient for systems analysis. Indeed, for a 
gaussian field these quantities are sufficient for a complete statistical 
characterization. 

Our most general results for a specific spectral-density function — 
namely, a power-law form — are given by Eqs. (III-33) through (III-37). 
They admit arbitrary incidence angles and a simple anisotropy, and no 
a-priori Fresnel approximation was made in deriving them. The results 
depend explicitly on the parameter p = c^k cos 0, where a is effectively 
the "outer" irregularity scale. 

When p is sufficiently large, a considerable simplification of the 
formulas can be realized. Note that p large is implied by the condition 
O' » \, which is expected to be valid even at HF frequencies. However, 
the effect of the large p terms depends on the spectral index y. We have 
found that when y > 1.5, p > 10 is generally sufficient to justify the 
Fresnel approximation (Section III-E). 
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When the Fresnel approximation is valid a general channel characteri- 

2 

sation is possible in terms of a single parameter Z = (\z sec 9)/2rav . 

That is, independent of all other parameters, the limiting behavior for 
large and small Z can be specified. For small Z, the phase-quad?’ature 
signal component is liiuch larger than the in-phase component. For lai’ge 
Z, the in-phase and phase -quadrature components are nearly equal. 

We have characterized these two behavior regions as the near and 
far zones respectively. The rate at which one proceeds from the near to 
the far zones as Z is increased from zero to infinity and the detailed 
behavior of the various quantities depend on the parameter values. 

For numerical computations we have used Eqs . (III-47) and (III-4S) 
with the parameters defined by Eqs. (III-49) through (III-52) . The 
quantities that are evaluated depend on the parameters listed in Table 
VI-1 below. The simplified formulas are valid whenever p is sufficiently 
large to justify the Fresnel approximation. 

The effects of the various parameters can be summarized as follows: 

• Ratio of Fresnel-zone area to scale parameter squared, Z. The 
value of Z determines to what extend free-space propagation 
has redistributed the scattered power from what it was in the 
vicinity of the scattering layer. When Z is sufficiently small 
the scattered power is nearly in phase quadrature with the 
coherent-signal component. When Z is sufficiently large the 
scattered power is nearly equally distributed between the 
in-phase and phase-quadrature components. 

• Layer thickness parameters, £,{{h/2y) contributes when G is 
finite]. The layer-thickness factor £ imposes a spatial filtex*- 
ing that suppresses spatial frequencies beyond x/£. Hence, its 
primary effect is for small Z, For example, it acts to keep 

finite when Z = 0 (Figure IV-4>. 

X 

• Two-frequency correlations /4f/f. The primary effect of finite 
Af/f is to make Z-dependent. Indeed, the functional 
dependence of % on Z and Af/f is nearly identical to that of 

with Z replaced by Z Af/f [Eq. (lV-48)]. It follows that 
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Tal3le VI-1 

PAEAJIETERS FOR MACHINE COMPUTATION 


Parameter 

Symbol 

Comments 

Scattering cross 
section 

2 a 2 

cr 1 sec Q 

2 

In Eq. ClV-47) is set 

equal to unity 
for convenience. 

Ratio of Fi’esnel zone area 
to scale parameter 
squai’ed 

Z 

A sec e 

2x0? 

Fractional frequency 
separation 

6f/f 

For t\vo-frequency 
correlations. Note 
that ]Af/±| < 1, 

Layer thiclmess 
parametei's 

r 

„ A Xh sec e 

" = ^2 
2 X 0 ! 

Ratio of layer thickness 
to scale parameter 

(L/ar) 

Only necessarj for 
oblique incidence 

Axial ratio and 
magnetic dip angle 

a, f 

Enters Eq. CIV-4S) as 
(a cos V + sin Y) 

Spectral index 

y 


Incidence angles 

9 , cp 

Theta is zenith angle. 
Phi is magnetic azimuth. 


the decorrelation of two frequencies is very slow for small 
Zj hence, the large cohei*ence handwidths that are characteris- 
tic of ionospheric scintillation. 


• Anisotropy, a, V, The anisotropy acts through p 4 (a cos V 

+ sin^ to reduce the rate of convergence to far-zone 

conditions as 0 increases from unity. Foi' isotropic irregulai*!- 
ties one observes the most rapid convergence. 

• Spectral index y. The spectral index acts effectively in the 
same way as the anisotropy. When y is large the spectrum has a 
sharp cutoff. Thus, the transition from near- to far-zone be- 
havior is abrupt. When y is small (near' unity), as it is 
generally observed to be, the transition is much slower. 
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• Oblique incidence, 9, tp. The dependences on incidence angles 
are generally quite complicated.. For example, the scattering 
cross section will exhibit a more complicated zenith-angle 
dependence than sec 0. The main impact is that for ionospheric 
modeling one can make large ex’rors if proper account for inci- 
dence angles is not taken. 

We have noted that for gaussinn fields, the complex covariances 
suffice to completely define the statistics . We have made use of this 
fact to calculate the intensity autocox*relation function. The results 
were applied to evaluate the wavelength dependence of the scintillation 
index and the two-frequency intensity correlations. We also outlined 
the method of computing phase statistics. 

The major limitation of our calculations stems from the restriction to 
weak scattering. We shall consider this limitation in detail below. However, 
we cannot expect a power-law spectral-density function with a constant 
spectral index to be accurate over an indefinite range. Indeed, the calcu- 
lated wavelength dependence of the scintillation index cannot account for the 
unexpectedly large scintillation observed at L-band and highei' frequencies . 

D. Consequences of the Weak-Scattering Restriction 

The most obvious problem with the first Born approximation is that 
it does not account for "extinction." That is, it gives no direct 
information about the coherent field component . We have- compensated for 
this deficiency in our calculations by reducing the incident field 
amplitude by the factox* Jl ~ cj^' so that enex'gy is conserved. 

Having done this , we observed no obvious breakdown of our intensity- 

statistics px’ocedui’es as the index became large (~0. G). However, 

our data base was limited. Nonetheless, the results do suggest that the 

/ 2 2 \ 

general neax'-zone stx’ucture of the field is likely to hold 

true for strong scattex’ing. At the present time we cannot be cex'tain 
of this. 
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We do know that we must ultimately accomit for the spectral hroaden- 
ing that accompanies multiple scattering and assess its effects as well 
as properly accounting for extinction. Unfortunately, there is no 
completely satisfactory method of proceeding. 

The parabolic-eqLuation approach combined with the Mai’kov approxi- 
mation is cm'rently gaining acceptance. The parabolic e<iUation is valid 
foi’ narrow-angle scattering, which is probably an acceptable I’estx’iction, 
However, the Markov approximation demands that the medlmn be "delta- 
functlon*'-corx’elated along the propagation direction. 

DeWolf (1972) has analyzed the full Box'n sei’ies sumiiiation with 
various approximations. His results show that the cox'X’elation terms that 
are neglected in the iMax*kov approximation can change the natux’e of the 
solution considex’ably. Moreover, anisotropy combined with the generally 
smaller ratio of the small-scale structure to wavelength for radio w-aves 
as Compared to optical waves, woi'ks to invalidate the Max’kov appx’oximation, 

We believe that the fact thai the fii-st-ox’der intensity statistics 
are most accui*ately reproduced by px'obability-density functions derived 
under the gauss Ian hypothesis is evidence that the Max'kov approximation 
is invalid for I’adio-wave scintillation. If it were valid, one would 
observe log-noimml statistics fox- weak scattex’ing. 

To relax the Markov approximation it appeal's that a dli’ect summation 
of the Box’n sei'ies is I'equired. One can achieve this indirectly by 
applying the weak-scattering results as an opei'atox* acting on an ai'bitrax’y 
incident field. The scattei'ing medium is divided into slabs that ax‘e 
large compared to the axial coi’X'elation distance. Then, the weak-scattex*- 
ing operator can be applied successively to each layei*. This is 
essentially the method described by Uscinski (1968). 

One would not expect the resuxLs to be analytically ti'actable. 

Howevei*, the method is amenable to machine computation. The necessary 
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integrations are similar to those that v/e have pei'ioi’med in deriving 
the numerical results presented in this i’eport. The disadvantages that 
result from not having an analytic solution are lax’gely offset by the 
fact that there are very few restrictions. 

B. Data Restrictions 

Scintillation data that nx’e analyzed for intensity statistics are 
limited by non-stationarity and noise contamination* Hence, for weakly 
scintillating signals one cannot obtain accurate statistics. The non- 
stationarlty limits the length of the data interval tliat can he i-epi’e- 
sented by a single histogram. Bor analysing higher-ordex’ statistics the 
problem is even more difficult. 

Diffraction effects complicate the intex’pretation of scintillation 
dataj as we have discussed in Section III-D. Howevex*, even when this is 
not the case, sample x’ates and noise limit our ability to analyze small- 
, scale structure. The ultimate breakdown of the assumption of statistical 
honogenelty Ximits our ability to accurately detei'mine the large-scale 
structure . At the present time the spatial wavenumber specti’Uin is not 
well defined below 1 km nox* above 50 km. 

The small-scale i-egion is probably important in explaining the 
scintillation that is being x*eported at L band and at G GBs. Thus, 
befox’e an accux'ate channel model can be constructed, a considex’able amount 
of careful data analysis will be required. 

To conclude , we x’eitei'ute our main x-esult — namely, that a gauss ian 
signal structure most accurately repx’oduces observed intensity histo- 
gi'atus. This fact to oui* knowlege has not previously been demonstx’ated . 
Aloreover, it is in conflict with the predictions of a lax’ge body of x'andom 
media propagation theory. 
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Appendix A 

THE SPECTRAL REPRESENTATION OP HOMOGENEOUS COMPLEX RANDOM PIELDS 
WITH APPLICATIONS TO PREE-SPACE PROPAGATION 
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Appendix A 

THE SPECTRAL EEPRESEETATIOIT OE EOMQGSNEOtlS COMPLEX! RANDOM FIELDS 
WITH APPLICATIONS TO FREE-SPACE PROPAGATION 

1. General Results 

In Section. Ill of this repo.rt the fundamental quantity is a scalar 
complex random field 

E(x,y,a) = Ej^Cx,y,z) + tE^Cx,y,z) . CA-1> 

For simplicity of notation we let E^ = ECx^,y^,z^) and = ECx" ,y'^ ,z" ) . 

Clearly, there are at least three spatial autocorrelation functions that 

characterize E— via. , {e^E'*^) , '(E E^) , and (E E^). Because of the layered 

’ R H I I R r 

model that we have used, E is homogeneous in any transverse plane. Hence, 
the correlations depend only on fix = x' - x'' , and fiy = y' - In 

general, the 2 ~dependence will involve both z* and z" . 

Consider Ej^. It admits a spectral representation of the form 

E Cx,y,z ) = rr exp{i<5C x + K y)3d5„CK ,K ja ) . (A-2) 

R o wU X y Rxyo 

Since is real, we must have 

dS.j3(K »(C J3 ) = dS*(ic ,k: ;z ) . (A-3) 

Rxyo Rxyo 


Now, has orthogonal CR ,h ) increments, 
R X y 

we can write 


In a purely formal notation, 


dtC dsc 
X y 


= fiCK' - K"n(K' - K'^)tp^(K ,R ja'.z'') ^ ^ 
R R X X y y X y 2x 2a 
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We exclude any discrete spectral co tiip clients . Because of Eq. CA-3), 
however, we can also wi'lte 

d5C _ die 

<c^S;dS") « 6CK^ + k")6CK' k")<o„CK ,IC ^ ♦ CA-S) 

R II X X y y X y 2« 2a 


We have a similar set of relations for E^ and d5j> as well ns for 
the cross-correlations (dS^dS^^*) and (clS^dS'l) . By comhining these 
results we can deduce the properties of the "transverse" specti’a! de- 
composition of E, 


ECx,y,j5 ) - ff exp[iCic x + w y)}«i?(ic ,w ;s ) 
0 lu ^ X y X y o 


CA-G) 


where 


dS = 


CA-7) 


We can write 


<dS^dS"*> = 6 Ck:' - - ic")«p„ U ,ic ]z ) 

X X y y x y o 


die die 

, X V 

2tr 2 it 


CA-8) 


dK , die 

<ds'ds"> = 6 Ck" -t- k")60c' t ()c ,k 

X X y y X y o 2st 2a 

ill 


where 


'I’r “ % 'fi *%i ■ 
E 


is the Fourier transform of R„ = (e'e"*), and 

E ^ * 


% = fE “ ®I '»Ie' 

E 


CA-9) 


is the Pouriei* transform of B^ = (E^e''), 
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The quantity tp^ is the ordinary angular spectral-density function. 

E . 

Indeed, since tp„_ = CD *, the angular spectral-density function is purely 
* RI ‘ IR 

■ 'O 

real. Moreover, since {E'*) S 0, it is also non-negative, although not 

necessarily symmetric. From the angular spectral-density alone, we cannot 

separate cpj^ and cp^. The additional quantity cp^ provides the necessary 

E 

information. As an example, we shall consider free-space propagation. 



By using Eq. (A-8) , we can immediately deduce the well-known result that 

R_C6x,6y,3) = R (6x,6yjz ) . CA-14) 

E E o 


However, cp^, and do depend on a. Indeed, by applying Eq. (A-9), 
we have 

tp„ ClC ,K ;z) = cp„ (k: ;z )expf2ik Cz - a )} . CA-15) 

X y ^B^ X y o z n ^ _ 

It follows that witr an obvious notation simplification. 
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I': 


f - J- 


J 




O 

tp^^Ca) s 9j^Ca^)cos [k^(a - a^)] 

2 * 

+ 9jCa^)sin "* 

Ca a )] CA-XB) 

' WJl O 5 o 


nud 


«p_(s) = {fl_(a )siii'^[k (a - 2 )] 

I o a 0 

+ tp„Ca )cos*^[k (a “• a )1 
0 a o 


- Ro{h) (a >}sinC2k Ca - a )] 
^ ivx 0 a o 


(A-XT) 


> “ C«Pi, ^ “ 'Pt ^ 3 s in[ 2k^ (s - a_^ ) 3 


»'H o no 


+ RG[tpj^jCa^)3«iost2k^Ca “ 




(A-XS ) 


Til the specitil case In which cpj^ = 0 at a = a , the resuXts 

reduce to 


2 

9j^Ca) s= tpj(a^)sin Ck^Ca “ 


(A-X9^ 


cpjCa) =: tp^(a^)cos‘"[k^ (a - 2 ^)] 


(A-20) 


tpp^Ca) = “ cp_(a >7 sln[ 2 k (a « z )] 
^Ri o 2 a o 


(A-21) 
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This special case was origihaliy derived by Bowhill (195Y). The results 
derived in Section III take this form when the layer thickness L is 
zero. This is the so-called ’'phase-changing" screen approximation, 
although "phase-duadrature" changing screen would be a more nearly 
correct description, A true phai5e=-changing screen has been analyzed by 
Itercier (1962 ) , 



101 



Appendix B 

THE TEMPORAL STRUCTURE OF THE SCATTERED FIELD 
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Appendix E 

THE TEMPORAL STRUCTURE OP THE SCATTERED FIELD 

The time variation of the scattered field arises from the relative 
motion of the medium, and the frame of reference as \vell as from temporal 
changes in the medium itself. Hence, in the measurement frame, 


iCK^Cx - u^t) 


+ K Cy - u t)]}d5 Ck a ;z) 

y y t X V 


y y ' t X y 

By appling Eqs , (A-8) and CA«0) \ve can derive the results 


(B-1) 


= ff exp{- i[K (6x - u 6t) 

X X 

dK dK 

+ K (6y ” u 6t)]}o (k ;6t) ^ ^ 


"z X y 


2k 2k 


. (B-2) 


In deriving Eq. (B-2) we assumed stationarity of the temporal 
variations , We introduce the temporal frequency variable 


2 t ^ = k u + k: u 

XX y y 


(B-3) 


Changing variables in Eq . CB-2) gives the result 
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Rl -A 


esp[ 21x176 1} 


'tH' 


2n7 - K u 

X X 

x’ u 


i6t j 


- f 


y / dK 


exp<2irl7 — }dj 


j f - K u \ u dK ( ■) 

tK-^s-y * ' 


CB-4) 


The form of Eq. (B-4) would be greatly simplified if $(k ,K ;6t) had the 

^ y 

form $(K^>Ky)n(6't) . It seems unlikely, however, that the decay mechanism 
is independent of spatial frequency or, equivalently, scale size. 

On the other hand, it is generally true that the lifetime of the 
irregularity structure is long compared to the transit time — that is, 
the irregularity structure is "frozen." Then, the integrals over K or 

X 

K^ are effectively 6t-independent , and they can be interpreted as the 

temporal spectral densities (power spectrum in the case of R ) of the 

E 

field fluctuations. 
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Appendix C 


POIVER-LAW SP!ECTRAL-DENSITy-FUNCTION FORMULAS FOR NORMAL INCIDEITCE 


In this appendix we consider evaluating the integral of Eq. (IV-40) 
for normal incidence with L - 0, By using hypergeometrio series, it is 
possible to obtain the following formulas for Re{B} and «sPm{B] : 


1 + Re[B} = D + C-D^^^C V^EA + B 1 


(C-1) 




(C-2) 


for Y and even Integer, and 


1 + Re{B} ='y^ S3 + C - B 1 (C-3) 

n y ^ 


-^mEB] ='^E + (^l/"^" C \^[A +B ] 


(C-4) 


for Y ail odd integer. 
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(- Y + 3/2) C- V -s- 5/2)... (- Y - 1/2 + 2n) 


A 

n 


2n + Y ~ 1 / 2 ^ 


C2n)I[3 


(2n + Y - 1 / 2 ) 



(C-7) 


B 

n 


, ^n. 2n + 

(-1) Z 


2n + Y 1/2 


(•^) 


(2n + Dip 


2n + Y + 1/2 


and 


C 

7 




1 ' 3 


• 5... (2Y - 3) 


(C-8) 


:c-9) 


where P Cx) is the O' Legendre polynoraical when cx is an integer or the 

O' 

appropriate Legendre function when ck' is non- integral. 


In formulas (C-7) and (C-S), the relations 




1 - r f 


i ! ‘ ! 

i 1 ■/ 

^.- 


I 



where K and E are standard elliptic integrals, can he used together with 

the I’ecuri’ence relations for P to evaluate the successive terms in the 

n 

series . 

The series in Eqs. (C-1), CC-2), (C-3), and (C-4) converge quite 
rapidly for small Z when y is not too large Csuy less than 6), and more 
slowly as Z is increased. Hence, they are most effective in the near 
zone. 
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Appendix D 

SPECTRAL-DENSITY FUNCTIONS FOR **S IMPLY" ANISOTROPIC MEDIA 
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Appendix D 

SPECTRAL-DENSITY FUNCTIONS FOR "SIMPLY" ANISOTROPIC MEDIA 


In general, the three-diraenstional spectral-density function is 


5Ch) = J(3)jR(6r) exp[ ik 


6 r} d6 r 


(D-1) 


where RCfir) is the spatial autocorrelation function. The z-axis is 

taken downward, with the x-z plane containing a unit vector u along 

B 

the earth's magnetic field. 


Consider a rotation about y into a new coordinate system (u ,u ,u ) 

12 3 

with u along ii , the elongation direction. In standard matrix notation 
o B 

we have 

u = M^6r CD-R") 

where 


/ cos ^ 

0 

-sin ijf\ 



Ml = ( 0 

1 

0 ) 

» 

<D-3) 

\ sin \if 

0 

cos ijf / 




The angle ijf is measured from z to u^. 

A second coordinate transformation to elongate the coordinate along 

the u direction produces a coordinate system with spherical symmetry 
o 

for the simple anisotropy we are considering — viz., elongation of the 
irregularities along a single direction. In the new system. 
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whei’e 


10 0 


M, = 0 I 0 


0 0b 


Now, changing coordinates in Eq. (D-1) gives 




•4 —1 — 1-t 

-./>=• ”2 V 

r)e dji 


whei’B a = b is the axial ratio. Let 


»+ 

S • ^ji = k • 61’ 


It follows that 


Hence , 


/ -1 -IN -4 

= (»>i » J s • 

= f M- H . 
12 


^1 ~ ^1 


Sn = K 
2 2 


Sg = + k^a sin + k^a cos 


,-»i2 *4 2 2 2 2 

jS| = S * S w k cos \{i 4 k sin ^ ijt sin ij; 

X O X <3 


2 2 2 2 2 2 2 
+ k ’■ sin + k a cos + 2k k a cos ijt sin i}t 
X o X o 


BCtlf) " ^^3] ■*■ \ BC^’O h 
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where 


and 


2 2 2 

BCiji) = sin ijj + a cos iji (D-12) 

C(ilj) - Cl - a®) cos i}r sin . CD-13) 


2 

We now note that $C|sj ) is a real symtneti’ic function. The quantity 
of interest Is 




= exp^ik 6zi 


.- 5lij /• ( 


2 2 a 2 

k BCijj ) + k -— + k 

3 1 B Cl) ) 2 


) if 


(D-14) 


The integral in Eq. CD-14) is purely real, as is the additional 
integral 


t3& 

/ $Ck ,k„;6a)d6z = §Ck •,k^,k ) Ik =0 . CD-15) 

«/ X 2 12 3*3 

— CO 


117 



J'BEGJIDIHG BAGE BLAtTS NOT FILMED 


Appendix E 

INTENSITY CORRELATIONS FOR COMPLEX GAUSSIAN FIELDS 


119 






■ftiMJiJDiiVG 


8LAi\k 




Appendix E 

INTENSITY CORRELATIONS FOR COMPLEX GAUSSIAN FIELDS 


In Section V we considered a general complex field of the form 
V = X + iY, where X and Y are jointly gaussian real fields. To make 
the deterministic components explicit, we write 


X = X + T , (E-l) 

Y = Y + T]y , (E-2) 

where X and Y have zero means. Note that we have not assumed statistical 
homogeneity. That is, T| and Tj can vary with position and/or time as 
well as the moments of X and Y. 

We let 


T 


R 4 <^0 , 

XX 

(E-3) 



R 4 . 

yy 

(E-4) 


and 

R k , 

xy 

CE-5) 



R 4 (X^Y> . 

yx . 

CE-6) 


When X = x' and Y = Y^ , we shall use the notation 


X 

<x^> 

(E-7) 

<11 

<Y^ 

(E-8) 

<11 

<XY> 

CE-9) 
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and 


2 2 2 
(j = a + cr 

X y 


(E-IO) 


Pinallyf we shall employ the auxiliary definitions 


R i = R + R + iCC - C ) 


XX yy xy yx 


(E-11) 


B ^ <VV> = R - R + i(C -f* C ) . 

XX yy xy yx 


(E-12) 


To proceed, we first note that 


<i) = <x^ + yS = 


CE-13) 


The quantitiy that we shall compute is tlie second central moment of 

9 

intensity (H^) - (1)“. By direct computation we have 

<II'> - <I>^ = 4[r + R -it + (R + R )T^ 1) 1 

L XX ’x yy y xy yx 'x 'yj 

+ + (Y^Y'^) -h <x\'^> + . CE-14) 


For gaussian processes it is easily shown that (see, for example, 
Papoulis, 1965, pp. 481-482), 

= (a^) + 2R^ CE-15) 

X XX 


. 2 ,2v 2 2 2 

<X Y^ > =; a cr + 2R 

X y xy 


CE-IG) 


with similar results for (Y^Y^% and (X^^Y*) . Hence, the intensity 
cori'elatlon can I 
moments T( _ and T) 


, 2 . 2 v 


cori'elatlon can he written in tei'ras of R , R , R , R , and the first 

XX yy xy yx 
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The final form can be simplified if we write B and (V) in polar form 


with 


and 


tan 2i* 4 


RefB} 


tan pf 4. — ^ 

'x 


(E-17) 


(E-IS) 


After some algebraic manipulations, we have the desired result, 


o of 


~Re[R} 

|fi| 1 

<II^> - (I> = 2a“(<l> -cr 

2 

+ cos 2Cc; “ sf) 


• 

- a 
1 |2- 

o'" 

4' 


^ Q 

4 

4 

* 


0 

O 



(E-19) 


^Ve note that both R and B are proportional to a . Hence, Eq. (E-19) 


is a quadi’atic equation in e * 
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